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jy  During  the- third  year  (FY'83)  a  great  deal  of _progress  has  been  made 
on  the  white-light  optical  information  processing  And  holography  research 
program.  this— period  ,  ~we— -have  completed  -the  work  of  a  broad  band  color 

image  deblurring  and  have— expanded  this  technique  to  the  restoration  of 

2-D  out-of-focused  color  photographic  images.  -We  have  also  developed  a^ 
new  technique  of  white-light  density  pseudocolor  encoder yfor  three  primary 
colors.  This  white-light  pseudocolor  encoder  is  very  cost,  effective  and  _ 
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[goffers  a  high  image  resolution,  which  would  be  an  excellent  alternation 
1  for  the  digital  counterpart.  We  have  also  in  this  pdriod  conducted  a 
measuring  technique  for  the  degree  of  coherence  in  the  Fourier  plane  of 
the  grating-based  white-light  signal  processor.  We  have  shown  that  high 
degree  of  coherence  is  achievable  with  a  signal  sampling  grating  it  the 
input  plane.  Thus  the  white-light  technique  is  capable  of  processing 
the  information  in  complex  amplitude  and  it  is  very  suitable  for  color 
signal  processing.  In  this  phase  of  research  we  have  also  studied  the 
effect  of  coherence  due  to  source  encoding,  signal  sampling  and  spectral 
band  filtering,  as  applied  to  the  white-light  signal  processing .'^^The 
effect  of  source  encoding  is  to  relax  the  constraints  of  a  physicalT"^ 
light  source,  and  the  effect  of  the  signal  sampling  is  to  improve  the 
degree  of  temporal  coherence  at  the  Fourier  plane.  Thus  the  signal  can 
be  carried  out  by  a  broad  spectral  filter.  In  this  period,  we  have 
investigated  the  noise  performance  of  a  white-light  signal  processor. 

The  noise  performance  under  the  temporally  incoherent  illumination  is 
quantitatively  analyzed.  We  have  shown  that  the  output  signal-to-noise 
ratio  improves  considerably  with  the  increasing  the  spectral  bandwidth 
of  the  light  source  employed.  In  addition,  we  have  also  evaluated  a 
dual-beam  encoding  technique  for  color  holographic  construction.  The 
color  hologram  image  can  be  easily  reproduced  by  a  white-light  processing 
technique.  This  encoding  technique  would  offer  the  simplicity  of  color 
hologram  generation  without  utilizing  a  reference  beam.  In  brief,  we 
have  shown  again  the  versatility  of  the  proposed  white-light  processing 
technique. 
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I.  Introduction 

In  the  past  year,  we  have  accomplished  several  major  research  works  on 
white-light  optical  information  processing  area.  Our  research  programs 
were  consistent  with  our  AFOSR  grant,  for  which  in  part  have  been  published 
in  various  refereed  journals  and  conference  proceedings,  and  have  been 
presented  to  several  scientific  conferences  (e.g.,  OSA,  SPIE).  Sample 
copies  of  these  articles  are  included  in  this  annual  report  in  the 
subsequent  sections,  to  provide  a  concise  documentation  of  our  research 
program.  In  the  following  sections,  we  shall  give  an  overview  of  our 
research  work  done  in  the  period  from  March  15,  1983  to  March  H,  1981!.  We 
will  detail  some  of  those  accomplished  works.  A  list  of  publications 
resulting  from  AFOSR' s  support  is  included  at  the  end  of  this  report. 

II.  Summary  and  Overview 

Although  optical  signal  processing  ha3  the  unique  abilities  of  (1) 
complex  amplitude  processing,  (2)  parallel  processing  capability,  (3)  high 
resolution,  (H)  speed,  and  (5)  color,  however  the  lack  of  versatility, 
stringent  processing  environment,  artifact  noise,  costly  light  sources, 
etc.,  those  are  the  major  impediments  for  the  slow  emerging  as  possible 
alternation  to  the  digital  counterpart.  Nevertheless,  there  are  techniques 
for  optical  signal  processing  that  can  successfully  compete  against  the 
electronic  and  digital  systems.  One  of  those  techniques  would  be  the 
white-light  signal  processing  technique.  The  advantages  of  the  white-light 
signal  processing  are:  (1)  It  eliminates  the  coherent  artifact,  (2)  the 
white-light  source  is  generally  inexpensive,  (3)  the  processing  environment 
is  rather  relaxed,  (*0  the  white-light  processor  is  easy  and  economical  to 
maintain,  and  (5)  the  processor  is  very  suitable  for  color  signal 
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We  have,  in  the  past  few  years,  developed  a  grating  based  white-light 
optical  signal  processing  technique  [1,2]  to  processing  the  optical  signal 
in  complex  amplitude.  We  have  illustrated  that  this  white-light  processing 
technique  is  capable  of  suppressing  the  coherent  artifact  as  an  incoherent 
processor  [3]  and  at  the  same  time  it  can  process  the  information  in 
complex  amplitude  as  a  coherent  processor  [4].  In  an  article  [5],  we  have 
shown,  that  smeared-photographic-image  can  be  deblurred  with  white-light 
processing  technique,  since  the  white-light  sources  possess  all  the  visible 
color  wavelengths,  it  is  particularly  suitable  for  color  image  processing: 
For  example,  the  color  image  deblurring  [6],  archival  storage  of  color  film 
[7],  pseudocolor  encoding  [8]  and  others.  We  have  also  developed  a  source 
encoding  concept  [9]  in  which  the  inability  of  an  extended  white-light 
source  may  be  alleviated.  That  is,  for  a  processing  operation,  a  source 
encoding  technique  may  be  generated  to  provide  a  specific  spatial  coherence 
function  that  is  suitable  for  the  processing,  e.g.,  the  image  subtraction 
with  encoded  incoherent  source  [10,11].  In  other  words,  a  broad  spatial 
coherence  function  may  not  be  needed  for  certain  optical  processing 
operations.  It  is  possible,  however,  to  obtain  a  reduced  spatial  coherence 
function  for  a  specific  processing  operation  by  spatially  encoding  the 
light  source.  We  have  also  extended  the  incoherent  image  subtraction 
technique  for  color  images  [11,12].  As  compared  with  the  results  obtained 
with  a  coherent  source,  we  see  that,  the  results  obtained  with  incoherent 
and  white-light  techniques  offer  a  higher  image  quality  than  those  obtained 
with  coherent  processsing  technique  [  1 3 ] - 

We  have,  also  in  the  past  two  years,  investigated  the  coherence 
requirement  for  this  white-light  optical  processing  technique  [14],  These 
basic  optical  processing  operations  are  evaluated.  We  have  shown  that  for 
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the  linear  Image  deblurring  problem,  the  spatial  coherence  requirement  is 
dependent  upon  the  source  size,  smeared  length,  and  spatial  frequency  of 
the  grating.  While  for  image  subtraction  problems,  the  coherence 
requirement  is  dependent  upon  the  ratio  of  slit  size  to  the  spading  of  the 
slits  of  a  source  encoding  mask,  the  slit  size,  the  spectral  bandwidth  of 
the  light  source  and  the  separation  of  the  input  object  transparencies. 

For  correlation  detection  [15],  the  requirement  for  temporal  coherence 
strongly  depends  on  the  spatial  frequency  and  the  spatial  extension  of  the 
target  (i.e.,  space  bandwidth  product).  However,  the  spatial  coherence 
requirement  depends  only  on  the  extension  of  the  target. 

We  have  also  formulated  an  apparent  transfer  function  of  this  proposed 
white-light  optical  processing  system  [16].  Although  these  formulas  show 
that  the  apparent  transfer  function  i3  dependent  upon  the  degree  of  spatial 
and  temporal  coherence,  there  is  actually  more  variability  in  the  spatial 
coherence.  We  note  that  the  obtained  formulas  may  also  be  used  as  a 
criterion  in  the  selection  of  source  size  and  spectral  bandwidth  of  an 
incoherent  light  source.  Thus  a  specific  optical  information  processing 
operation  can  be  carried  out  with  an  incoherent  source. 

We  have  also  in  this  period  evaluated  the  (white-light  transmission) 
rainbow  holographic  aberrations  and  its  bandwidth  requirements  [17].  We 
have  shown  the  conditions  for  the  elimination  of  the  five  primary  rainbow 
holographic  aberrations.  These  conditions  may  be  useful  for  the 
application  of  obtaining  a  high  quality  rainbow  hologram  image.  In  term3 
of  bandwidth  requirement,  we  have  shown  that  the  bandwidth  requirement  for 
a  rainbow  holographic  construction  is  usually  several  orders  lower  than 
that  of  a  conventional  holographic  process,  hence  a  lower  resolution 
recording  medium  can  generally  be  used. 
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We  have  also  conducted  a  research  on  the  detection  of  phase  object 
variation  through  color  encoding  with  encoded  extended  incoherent  sources 
[183.  This  technique  provides  fine  detail  of  the  object  phase  variation 
between  fringes,  including  both  the  positive  and  negative  phase  variations. 

In  the  following  sections,  we  shall  highlight  some  of  the  research 
done  during  the  past  year  effort  on  white-light  signal  processing. 

2.1  Broad  Spectral  Band  Color  Image  Deblurring  (Section  III) 

In  the  past  twelve  months,  we  have  conducted  a  broad  spectral  band 
color  image  deblurring  [193.  The  technique  utilizes  a  fan-shape  broad 
spectral  band  deblurring  filter  to  compensate  the  scale  variation  of  the 
smeared  Fourier  spectra.  This  technique  i3  particularly  suitable  for 
linear  smeared  color  image  deblurring. 

2.2  Noise  Performance 

In  this  period,  we  have  also  quantitatively  analyzed  the  noise 
performance  for  a  white-light  signal  processor.  The  analysis  is  broken 
down  into  two  major  parts;  the  effects  due  to  temporal  coherence  and  due  to 
spatial  coherence.  We  have  completed  the  noise  performance  due  to  the 
temporally  partially  coherent  illumination,  for  which  a  paper  is  submitted 
for  publication  [203. 

2.3  Pseudocolor  Encoding  with  Three  Primary  Colors  (Section  IV). 

We  have  also  in  this  period  developed  a  white-light  density 

pseudocolor  encoder  for  three  primary  colors  [213.  The  advantages  of  this 
technique  are;  it  is  very  co3t  effective  arid  offers  a  high  image 
resolution,  as  compared  with  the  digital  technique. 

2.4  Partial  Coherent  Measurement  (Section  V). 

We  have,  in  this  phase,  accomplished  a  coherent  measurement  for  our 
white-light  optical  signal  processor  [ 22 3.  The  results  show  that  the 


degree  of  coherence  in  this  Fourier  plane  increases  as  the  spatial 
frequency  of  the  sampling  grating  increases.  However,  the  improvement  in 
coherence  is  somewhat  more  effective  in  the  direction  perpendicular  to 
light  dispersion.  Since  the  white-light  processor  is  capable  of  processing 
complex  signal  with  entire  spectral  band  of  the  light  source,  it  is 
suitable  for  color  signal  processing. 

2.5  Restoration  of  Out-of-Focused  Color  Image  (Section  VI). 

In  this  period,  we  have  also  extended  the  color  image  deblurring  to 
2-D  out-of-focused  color  photographic  images  [23].  Final  results  of  this 
restoration  technique  have  been  obtained. 

2.6  Source  Encoding,  Sampling  and  Spectral  Band  Filtering  (Section  VII). 

We  have  also  in  the  past  year  developed  a  general  concept  on  source 

encoding,  signal  sampling  and  spectral  band  filtering  for  a  white-light 
processor  [24]..  A  similar  paper  was  also  presented  to  the  10th 
International  Optical  Computing  Conference  at  M.I.T.  on  April  6-8  [25]. 

2.7  Advances  in  White-Light  Optical  Signal  Processing  (Section  VIII). 

We  have  also  in  this  period  reported  a  paper  on  the  recent  advances  in 
white-light  to  the  proceedings  on  Optical  Information  Processing  Conference 
II,  NASA  Langley  Research  Center  [26]. 

2.8  Remarks 

We  have  also,  in  this  period,  investigated  the  real-time  and  computer 
controllable  processing  capability  as  applied  to  white-light  signal 
processing.  We  are  also  currently  developing  a  computer  generated  filter 
algorithm  for  our  white-light  signal  processor.  Interesting  results  would 
be  surfaced  in  the  near  future.  In  short,  the  white-light  signal 
processing  research  program,  supported  by  AFOSR ,  is  conducted  extremely 
well  as  proposed.  As  it  can  be  seen,  several  significant  results  have  been 
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documented  in  open  literatures  available  for  the  interested  technical  and 
research  staffs,  for  example,  as  listed  in  Section  III  to  VIII. 

2.9  Future  Research 

Our  aim  in  the  following  years  is  toward  the  following  goal's: 

1.  The  synthesize  a  low  cost  computer  controllable  white-light 
pseudocolor  encoder. 

2.  To  carry  out  tho  white-light  programmable  pseudocolor  encoder  for 
multi spectral  band  satellite  pictures,  x-ray  transparencies,  scanning 
electron  micrograph,  etc. 

3.  Measurement  of  the  noise  performance  of  the  white-light  signal 
processor. 

i*.  To  develop  a  real-time  programmable  processing  capability  for  the 
white-light  signal  processing  technique. 

5.  To  carry  out  various  real-time  signal  processing  applications. 

6.  To  develop  a  real-time  color  signal  processing  capability. 

7.  Provide  experimental  demonstrations  and  applications  of  the 
principles  and  processing  operation  that  we  proposed. 

8.  To  develop  techniques  of  synthesizing  a  broad  spectral  band 
complex  matched  filter  for  the  white-light  processing  system. 

9.  To  develop  a  computer  generated  spatial  filters  program  that  is 
suitable  for  broad  spectral  band  white-light  processing. 

10.  To  develop  a  technique  of  utilizing  real-time  magneto-optics 
spatial  light  modulator  for  filter  synthesis,  as  applied  to  white-light. 

11.  Improve  the  resolution,  system  performance,  and  signal- to-noise 
ratio  of  the  proposed  system. 

12.  To  develop  a  generalized  source  encoding  technique  for  the 
proposed  white-light  processor. 
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13.  To  develop  a  sun-light  processing  capability. 

111.  Evaluate  the  merits,  constraints,  and  limitations  of  the  proposed 
system. 

15.  To  develop  techniques  of  generating  and  reconstructing  hologram 
philosophies  of  the  US  Air  Force  and  DOD  needs. 

The  modern  tactical  and  strategic  Air  Force  engagement  scenarios 
assume  rapid,  effective,  secure,  independent  processing  and  communication 
among  many  ground  and  airborne  stations.  We  believe  that  the  proposed  work 
will  have  a  profound  and  direct  effort  on  Air  Force  processing  system  as 
well  as  provide  substantial  "fall  out"  benefit  in  such  areas  as  radar, 
sonar,  tracking,  targeting,  surveillance,  inspection,  and  many  others. 
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Broad  spectral  band  color  Image  deblurring 


T.  H.  Chao,  S.  L.  Zhuang,  S.  Z.  Mao,  and  F.  T.  S.  Yu 


A  broadband  while- light  processing  technique  for  smeared  color  photographic  image  dchlurring  isdescrilied. 
The  technique  utilizes  a  diffraction  grating  method  to  disperse  the  smeared  image  spectra  in  the  Fourier 
plane  so  that  the  entire  s|tcetrnl  band  of  the  while-light  source  can  l>c  utilized  lor  the  dchlurring.  In  Ibis 
paper  tile  technique  of  synthesizing  a  fnn-shn|M‘  ly|M-  complex  dchlurring  filler  In  accommodate  wavelength 
variation  is  presented.  Kx|M-rimenlal  results  showed  that  this  broad  spectral  baud  processing  technique  of¬ 
fers  an  excellent  coherent  artifact  noise  suppression,  and  the  technique  is  particularly  suitable  for  color 
image  deblurring.  Experimental  demonstrations  and  comparisons  with  the  narrowband  and  coherent  de¬ 
blurring  are  nlso  provided. 


I.  Introduction 

Restoration  of  smeared  photographic  images  lifts  lout; 
been  an  interesting  and  important  application  in  optical 
processing.'  7  In  image  dchlurring  much  of  the  effort 
has  lieen  devoted  to  applying  inverse  filtering  concepts 
to  the  image  restoration.  As  those  works  evolved,  two 
problems  are  still  of  intense  interest;  namely,  the  co¬ 
herent  artifact  reduction  and  color  image  dehlurring. 
Although  Wiener  filters  had  been  applied  in  coherent 
processors  hv  several  invest  ig.at  ions'  11  for  noise  reduc¬ 
tion,  they  do  not  suppress  the  inherent  coherent  artifact 
noise  in  the  processing  system.  Recently,  Yang  and 
Leith7  proposed  a  spatial  domain  deconvolution  tech¬ 
nique  for  image  dehlurring.  They  used  an  extended 
incoherent  line  source  for  dehlurring,  and  the  coherent 
noise  was  remarkably  reduced.  However,  their  tech¬ 
nique  is  only  suitable  for  processing  monochrome 
blurred  images.  In  previous  papers8  '0  we  presented 
a  white-light  processing  technique  for  linearly  smeared 
image  dehlurring.  We  have  shown  that  the  white-light 
image  dehlurring  technique  is  capable  of  eliminating  the 
coherent  artifact  noise  and  is  suitable  for  color  image 
dehlurring.  However,  the  results  that  we  obtained  were 
primarily  restricted  to  the  narrow  spectral  band  de¬ 
hlurring  concept. 

In  this  paper  we  shall  extend  the  image  dchlurring 
technique  to  the  entire  broad  spectral  hand  of  the 


white-light  source.  To  obtain  the  broadband  de¬ 
hlurring  effect ,  a  fan-shaped  spat  ial  filter  to  compensate 
the  scale  of  the  Fourier  spectra  should  he  ut  ilized  at  the 
Fourier  plane.  'Pin*  technique  of  synthesizing  the 
fan-shaped  broadband  dehlurring  filter  is  given.  The 
coherence  requirements  for  the  while-light  image  dc¬ 
hlurring  are  illustrated.  Experimental  demonstrations 
with  the  comparison  of  narrowband  dehlurring  and 
coherent  technique  are  provided. 

II.  Broadband  Image  Deblurring 

We  shall  now  discuss  a  broadband  image  dehlurring 
technique  utilizing  the  entire  spectral  hand  of  a  white- 
light  source.  Let  this  linear  smeared  imag'  he  given, 


*(x,v)  «  *(x,.v)  •  reel  •  ’  (1) 

where  s(x,y)  and  .s(x,y)  are  the  smeared  and  unsmeared 
images, 


-\  * 

W)  0. 
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othprwise. 
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and  W  is  the  smeared  length. 

Let  us  insert  the  smeared  image  transparency  of  F,q. 
( I )  into  the  input  plane  /’  i  of  a  white  light  optical  pro¬ 
cessor  as  shown  in  Fig.  1 .  The  complex  light  distribu¬ 
tion  for  every  wavelength  X  at  the  hack  focal  length  of 
the  transform  lens  would  be 

/i(of,0;X)  ■=  C  Jf  s(x,v,A)  exp(ipox) 


X  exp  |-  i  —  (xe  +  ydlj  dxdv.  (3) 

where  p<i  is  the  angular  spatial  frequency  of  the  phase 
grating,  or  =  \(Xf)/2v\p  and  fi  =  |(X/)/2ir|t/  represent  the 
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where 


h|x*I ini  coordinate  system  of  Fourier  plane  / *-j,  (/>,</)  is 
the  corres|sinding  angular  spatial  frequency  coordinate 
system.  /  is  the  ftsiil  length  of  I  lie  achromatic  transform 
Ions,  mul  ( ’  is  an  iipproprialo  complex  constant.  Thus, 
Kq.  (3)  can  be  written  as 

(a-^po^).  X4) 

S|a  -  (X//2*)pn,0|  -  S  jrr  -  pp,/jJ  sine  flj  (5) 

is  I  lie  linear  smeared  image  spect  rum. 

Since  scale  of  the  signal  spectrum  is  proportional  to 
the  wavelength  of  the  source,  the  corresponding 
signal  spectra  would  smear  into  a  fan-shaped  rainlmw 
color  ns  ran  lie  seen  from  K<|.  (ft).  In  ol her  words,  the 
top  (i.e.,  the  wider  region)  of  the  smeared  spectra  is  in 
red  and  the  bottom  is  in  violet. 

Let  us  assume  that  a  fan-shaped  broad  spectral  band 
deblurring  filter  (i.e.,  a  broadband  inverse  filter)  to  ac¬ 
commodate  the  variation  of  the  scale  of  the  signal 
spectra  is  available.  This  fan-shaped  filter  is  described 
in  the  following  equation: 

HU, fiM  -  X  (,.  -  jVd)  (J  («•<•«  (J)  «P  (> «  ~  *)|  4'  > 


(6) 


In  image  deblurring  we  would  insert  this  dehlurring 
filler  of  Eq.  (ft)  in  the  spatial  frequency  plane  of  P?. 
The  complex  light  distribution  for  every  A  at  the  output 
image  plane  P. t  can  he  written  as 


g(.r „v:A)  *  F  1  sy'  ~  //,n.d;X) 


(7) 


Fig.  I.  White-light  processor  for  smeared  color  image  dehlurring: 
7.  while-light  point  source;  .^(*,v),  smeared  color  image  transparency; 
T(rl.  diffraction  grating;  !.\  and  l.-t,  achromatic  transform  lenses; 
lfln,il;X|.  hroad  spectral  hand  dehlurring  filler. 
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Fig.  2.  Plots  of  l  lm  dehlurring  widlh  All’  ns  a  function  of  the  s|NS'trnl 
Isuulwidlh  of  Us*  light  source  AX  lor  variisis  vnliHst  of  nmennsl  length 

IV. 

where  F~l  denotes  the  inverse  Fourier  transform;  by 
substituting  Kqs.  (ft)  and  (f>)  into  Kq.  (7),  wc  have 


/?(x.y;X)  =  *tx,y)  expfip.rf). 


(8) 


which  is  independent  of  the  wavelength  of  the  light 
source.  The  resultant  output  intensity  distribution  can 
be  shown  as 


7(x,y)  ’ 


f  |g(xA,;X)|2dX 

•/ AX 


AX|*(x^)|2 


(9) 


which  is  proportional  to  the  entire  spectral  bandwidth 
AX  of  the  white-light  source.  Thus  we  see  that  this 
proposed  white-light  dehlurring  technique  is  capable 
of  processing  the  information  with  the  entire  visible 
spectral  band,  and  it  is  very  suitable  for  the  application 
to  color-image  dehlurring.  Since  the  integration  of  Eq. 
(9)  is  taken  from  the  entire  spectral  band  of  the  white- 
light  source,  the  coherent  artifact  noise  in  principle  can 
be  eliminated. 

III.  Coherence  Requirement 

Although  this  proposed  dehlurring  technique  11tili7.es 
n  while-light  source,  the  processing  is  operated  in  a 
partially  coherent  mode.  It  is,  therefore,  our  aim  in  this 
section  to  discuss  the  basic  coherence  requirement  for 
this  proposed  color  image  deblurring  technique. 

In  a  previous  paper*1  we  obtained  the  coherence  re- 
quirerr  ents  for  a  partially  coherent  optical  processor. 
Several  of  those  fruitful  results  can  be  applied  to  our 
proposed  white-light  image  dehlurring  system. 

We  shall  first  discuss  the  temporal  coherence  re¬ 
quirement  for  the  image  dehlurring.  We  shall  use  the 
results  obtained  in  our  previous  article  as  shown  in  Fig. 
2.  This  figure  shows  the  plots  of  the  spectral  width  AX 
requirement  of  the  light  source  (i.e.,  equivalent  to  t  he 
narrow  spectral  width  of  the  dehlurring  filter)  ns  a 
function  of  dclilurred  length  (A IV)  for  various  values 
of  smeared  length  W.  The  wavelength  spread  across 
a  narrow  spectral  band  filter  centered  at  wavelength 
A,,™  is 


AX  *  X|,|4 A/>/j)«|, /)»  »  A p. 


U0) 


or 


til) 


Ap  AX 
Ihi  -IX,, 

where  A p  is  the  angular  spatial  frequency  limit  of  the 
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Fig.  3.  Plots  of  the  deblurring  width  as  a  function  of  the  source  size 
for  various  values  of  the  smeared  length  IV. 


Table  II.  Effect  of  Spatial  Coherence  Requirement 


N.  AW 

\  .\h 

\us  mnr\ 

w\  \ 
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0.(1 

0.92 

1  mm 

0.1 

0.1 8 

0.20 

0.40 
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0.05 

0.08 

0.12 

0.18 

input  blurred  object  transparency,  and  />u  is  the  angular 
spat  ial  frequency  of  t  lie  phase  grid  inj*.  The  tcni|mrnl 
coherence  requirement  for  a  narrow  spectral  hand  de¬ 
hlurring  filter  for  wavelength  5461  A  is  tabulated  in 
Table  I. 


We  shall  now  discuss  the  spatial  coherence  require-  - 
men  I .  Since  the  dehlurring  operation  nets  on  the 
smeared  length,  the  coherence  requirement  is  depen¬ 
dent  on  the  width  of  the  light  source.  Let  us  now  use 
the  plots  of  dehlnrring  width  as  a  function  of  the  source 
size  as  shown  in  Fig.  3.  From  this  figure  we  see  that  the 
dehlnrring  effect  (i.e.,  AW)  ceases  when  the  source 
width  AS  reaches  a  critical  width  AS.  We  shall  now 
summarize  the  results  of  the  spatial  coherence  re¬ 
quirement  in  Table  II.  From  this  table  we  see  that  the 
higher  the  degree  of  dehlnrring  (i.e.,  smnller  AW),  the 
more  critical  the  spatial  coherence  (i.e.,  smaller  the  AS) 
is  required.  It  is,  therefore,  one  of  the  prices  we  paid 
for  a  higher  degree  of  dehlurring. 

IV.  Broadband  Deblurrlng  Filter  Synthesis 

We  shall  now  briefly  describe  (he  synthesis  of  a  fan¬ 
shaped  (i.e.,  broad  spectral  hand)  dehlurring  filter.  The 
synthesis  is  a  combination  of  an  absorptive-amplitude 
filter  and  a  phase  filter.  A  fan-shaped  phase  filter  is 
composed  of  several  slanted  bar-type  phase  objects  as 
illustrated  in  Fig.  4.  Kach  phase  bar  would  give  rise  to 
specific  ir  phase  retardation  fora  predescribed  disper¬ 
sion  of  rainbow  color  wavelength.  We  note  that  the 
height  of  dehlnrring  filter  is,  of  course,  dependent  on  the 
grat  ing  frequency  pn  at  the  input  plane.  The  period¬ 
icity  of  the  deblurred  filter  is  certainly  determined  by 
the  smeared  length  of  the  blurred  object,  and  the  width 
of  the  filter  defines  the  degree  of  dehlurring.1'2  In 
constructing  a  broad  spectral  hand  phase  deblurring 
filter,  we  utilize  a  vacuum  deposition  technique.  It  can 
he  accomplished  by  dcftositing  t  he  magnesium  fluoride 
(MgFj)  on  this  surface  of  an  optical  Hal  glass  stilts!  rate. 
In  this  technique,  a  blocking  mask  of  a  fan -shaped  bar 
pattern  as  shown  in  Fig.  4  is  used  for  the  vacuum  de¬ 
position.  The  MgF^  vapor  is  deposited  through  this 
blocking  mask,  together  with  a  linear  moving  covering 


Although  the  temporal  requirement  is  based  on  a  plate,  from  top  to  bottom  as  illustrated  in  Fig.  4. 
narrow  spectral  hand  analysis,  it  can  be  extended  to  a  The  thickness  of  the  deposited  coating  can  be  deter- 
broad  spectral  hand  operation.  For  example,  a  broad  mined  by  the  following  equation: 

spectral  band  filter  (e.g.,  fan-shaped  deblurring  filter)  ^ 

can  be  considered  as  a  summation  of  a  sequence  of  ~ _  ~  •  (12) 

narrowband  spatial  filters  of  various  wavelengths.  The 
deblurred  image  is  the  result  of  the  superposition  of  the 


mutually  incoherent  light  fields  derived  from  the  nar¬ 
rowband  filters.  Thus,  the  temporal  coherence  re¬ 
quirement  shown  in  Table  I  can  also  be  applied  to  a 
broad  spectral  hand  filtering.  As  an  example,  if 
(A /))//>(,  =  0.012and  the  spectral  bandwidth  of  the  light 
source  is  3000  A,  the  broad  spectral  hand  dehlurring 
filter  is  approximately  equal  to  the  sum  of  eleven  narrow 
spectral  hand  filters.  The  mean  spectral  bandwidth  A\ 
of  those  filters  is  270  A,  and  the  dehlurring  ratio 
(A  W)/W  is  1/20.  In  other  words,  it  is  possible  tosyn- 
thesi/.e  a  fan-shaped  type  spatial  filter  to  compensate 
with  the  scale  variation  of  the  smeared  Fourier  spectra 
in  the  spatial  frequency  plane  so  that  the  dehlurring 
takes  place  with  the  entire  spectral  band  of  the  white- 
light  source.  Since  the  broadband  deblurring  utilizes 
the  whole  visible  spectrum  of  the  light  source,  it  is 
particularly  suitable  for  color  image  deblurring. 
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Fig.  4.  Phase  filter  mask  for  MrF2  vapor  deposition. 
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where  n  in  the  ref  rad ive  index  iift.hu coaling  material. 
This  coating  thickness  is  linearly  proportional  to  the 
dispersion  of  the  illuminating  wavelength.  We  note 
(hat  a  strictly  linear  control  of  coaling  thicknesses  is 
very  essential.  The  advantage  of  this  phase  ty|>e  filter 
is  to  improve  the  transmission  efficiency,  since  the 
overall  deblurring  filter  is,  in  general,  highly  absorp¬ 
tive. 

In  principle,  it  is  a  straightforward  method  to  syn¬ 
thesize  a  fan-shaped  amplitude  filter.  The  synthesis 
can  be  accomplished  by  inserting  a  slit  aperture  of  a  slit 
width  ecpinl  to  the  smeared  length  of  the  blurred  image 
at  the  input  plane /*»  of  the  white-light  optical  processor 
shown  in  Kig.  1.  The  size  of  the  white-light  source 
should  lie  adequately  small  under  I  lie  spatial  coherence 
regime  to  obtain  a  smeared  sine  factor  (i.e.,  smeared 
Fourier  spectra  of  the  slit  aperture)  in  the  Fourier  plane. 
An  amplitude  filter  can  then  he  synthesized  by  simply 
recording  this  smeared  sine  factor  on  a  photographic 
plate.  If  the  film-gamma  of  the  recorded  plate  is  con¬ 
trolled  to  alsiut  unity  (i.e.,  y  -  1),  the  amplitude 
transmittance  or  the  recorded  plate  is  equivalent  to  that 
of  the  desired  fan-shaped  amplitude  filter.  However, 
in  practice,  a  fan  shaped  type  amplitude  filler  is  not 
that  easy  to  synthesize  due  to  three  primary  reasons. 
First,  if  a  very  small  source  size  is  required  for  the  filler 
synthesis,  it  usually  takes  a  longer  exposure  time,  for 
example,  if  Kodak  64DF  plate  (a  low-speed  Him)  is  used. 
Second,  the  spectral  response  of  the  recording  plate  is 
generally  not  uniform  for  all  visible  wavelengths.  The 
recorded  filter  would  produce  uneven  transmittance  in 
the  direction  of  the  smeared  color  spectra.  The  effect 
of  the  transmittance  variation  of  the  filter  would  affect 
the  fidelity  of  color  reproduction  and  the  degree  of 
restoration.  Third.it  is  difficult  to  synthesize  a  side¬ 
band  amplitude  tiller,  since  the  dynamic  range  of  the 
photographic  film  is  very  limited. 

There  is  an  alternative  technique  of  generating  a 
fan-shaped  amplitude  filter  with  coherent  illumination 
ns  shown  in  Fig.  5.  The  purpose  of  using  a  curved-slit 
aperture  is  to  accommodate  the  scab;  variation  of  the 
amplitude  filler.  The  expression  of  the  curved-slit 
aperture  can  lie  written 

rflr.v.Xl  ■  ml  (  — - — |  f)  (r - -  .  ( I M ) 

\W  A„/Xj  1  2*  / 

where  A(1  is  the  wavelength  of  the  coherence  source,  W 
is  the  smeared  length  of  the  hlurred  image,  f  is  the  focal 
length  of  the  cylindrical  transform  lens,  and  \  is  the 
wavelength  of  the  white-light  source. 

The  corresponding  Fourier  I ranslormal ion  of  the 
curved -slit  aperture  can  he  shown  as 

/J(i».d;A)  *  sine  •  x|i«  -  ’  04) 

where  *  denotes  the  convolution  operation. 

It  is  clear  now  that  a  photographic  recording  of  the 
spectra  shown  in  Eq.  (14)  would  produce  a  desirable 
fan-shaped  amplitude  filter  for  deblurring.  In  synthesis 
of  this  broadband  amplitude  filter,  a  He-Ne  laser,  with 
a  rotating  ground  glass  to  reduce  the  artifact  noise,  is 


-  i 


Fig.  ft.  ( of  iiinplihiOc  Hllrr  .S'  mm i<m  Itrnttinlir  |»lntt<»  wnvr; 

(7#,  rylitiilrintl  Irmmlnrm  Irns. 

used  ns  n  coherent,  source.  Kodak  Ft  I  plat  e  is  used  for 
the  recording  plaLe,  and  a  G-min  developing  lime  in  a 
I’OTA  developer  at  24 •(■  is  used  to  control  the  film 
gamma  to  alioul  unity.  The  spectral  wavelength  limits 
are  rhosen  from  4000  to  7000  A.  Within  the  dynamic 
range  of  the  recording  film,  five  sidelobes  of  ~.'I00  1 
dynamic  range  (or  n  density  range  of  2..r>)  are  recorded 
with  good  accuracy.  The  fail  shaped  amplitude  filter 
obtained  is  tested  with  satisfactory  results. 

V.  Experimental  Results 

In  this  section  we  shall  provide  a  few  experimental 
results  of  image  dehlurring  utilizing  a  broadband 
white-light  source.  In  our  experiments,  a  75-W  xenon 
arc  lamp  with  a  200-pm  pinhole  is  used  as  a  broadband 
white-light  source.  A  phase  grating  of  I  MO  lines/mm 
with  2.V’}.  diffraction  efficiency  at  each  first  order  dif¬ 
fraction  is  used  at  the  input  plane.  An//H  transform 
lens  wit  h  M00  mm  focal  lengl  h  is  used  for  image  Fourier 
transformation. 

We  shall  first  demonstrate  the  effect  of  the  broad¬ 
band  deblurring  ns  compared  with  the  narrowband  and 
the  result  obtained  with  coherent  source.  For  sim¬ 
plicity  of  illustrations,  we  use  a  set  of  linear  blurred  al¬ 
phabets  as  input  objects  as  shown  in  Fig.  6(a).  The 
smeared  length  is  ~0.,ri  nun.  Figure  6(b)  shows  the 
deblurred  image  obtained  with  Ibis  broadband  de¬ 
blurring  technique, and  lhes|ieclrnl  bandwidth  is~M000 
A  under  white-light  illumination.  Figure  6(c)  is  the 
result  obtained  with  a  narrow  spectral  band  deblurring 
filter  of  ~AX  =  100  A  centered  at  0M2H  A.  Figure  6(d) 
is  the  deblurred  image  obtained  wit  h  a  I  le  Ne  coherent 
source.  In  the  comparison  of  these  results,  we  see  that 
the  results  obtained  with  a  broad  spectral  band  white- 
light  source  offers  a  higher  deblurred  image  quality;  for 
example,  the  coherent  artifact  noise  is  substantially 
suppressed,  and  the  deblurred  image  appears  to  be 
sharper  than  the  one  obtained  with  a  narrowband 
case. 

We  shall  now  experimentally  demonstrate  the  ca¬ 
pability  of  the  white-light  technique  for  color  images. 
Figure  7(a)  shows  a  black-and-white  color  blurred  image 
of  a  building  due  to  linear  motion  as  an  input  color 
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PENN  PENN 
STATE  STATE 

(a)  «» 


PENN  PENN 

STATE  STATE 

(c)  (d) 

Fig.  fi.  Smeared  imago  rcslornlion  of  (he words  I’KNN STATK:  (n) 
smeared  linage.  (Ii)  deltlurred  image  obtained  with  broadband 
white-light  sourer,  (r)  debhirred  imago  obtained  wilh  narrow  H|M>ctral 
band  while-light  source,  and  (d)  dcblurred  image  obtained  with  co¬ 
herent  source. 


(b) 

Fig.  7.  Continuous-tone  color  image  deblurring:  (a)  a  black-and- 
white  picture  of  a  smeared  color  photograph  of  a  building,  (b)  a 
black-and-white  picture  of  the  deblurred  color  image. 


object.  From  this  figure  we  see  that  the  white  window 
frames,  the  front  doors,  bushes,  two  white  Itemns,  trees, 
elf.  are  severely  smeared.  Figure  7(b)  shows  the  de- 
blurring  result  that  we  have  obtained  with  this  white- 
light  processing  technique.  From  this  figure,  we  have 
seen  that  the  color  reproduction  is  rather  faithful  and 
the  deblurred  effect  is  sjtectaculariy  good;  for  example, 
the  window  frames,  the  bushes,  the  beams,  the  front 
doors,  the  trees,  etc.  can  be  clearly  identified. 

We  would  now  provide  another  more  striking  exam¬ 
ple  of  the  color  image  deblurring  wilh  Ibis  while-light 
processing  technique.  Figure  8(a)  shows  a  black-and- 
white  linear-motion  blurred  picture  of  an  F-16  fighter 
plane.  The  body  of  this  fighter  plane  is  painted  in 
bluc-niid-whilc  colors,  (lie  wings  arc  mostly  painted  red, 
the  tail  is  hlue-and-whitc,  and  the  ground  terrain  is 
generally  a  bluish  color.  From  (his  figure  we  see  that 
the  letters  on  the  body  on  one  of  the  wings  and  on  this 
side  of  the  tail  are  smeared  beyond  recognition.  The 
details  of  the  missiles  at  the  tips  of  the  wings  are  lost. 


(b) 

Fig.  8.  Color  image  deblurring:  (a)  a  black-and-white  picture  of  a 
smeared  color  image  of  an  F-16  fighter  plane,  (b)  a  black-and  white 
picture  of  the  deblurred  color  image. 
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The  star  symbols  at  the  tail  end  of  the  body  and  on  the 
top  of  the  winus  are  badly  distorted.  The  features  of 
ground  terrain  are  obscured.  Figure  8(b)  shows  the 
color  image  deblurring  result  that  we  have  obtained 
from  Fig.  8(a)  with  the  proposed  white-light  deblurring 
technique.  From  this  deblurred  result,  the  letters 
USAF  on  the  wing  and  YF-16  on  the  side  of  the  tail  can 
be  clearly  seen.  The  words  U.S.  AIR  FORCE  may  be 
recognized.  The  star  symbol  on  the  wing  can  be  clearly 
identified;  however,  the  one  on  the  body  is  rather  ob¬ 
scured.  Undoubtedly,  the  missiles  at  the  tips  of  the 
wings  can  be  seen,  and  the  pilot  in  the  cockpit  is  quite 
visible.  The  overall  shape  of  the  ent  ire  airplane  is  more 
distinctive  than  the  blurred  one.  Moreover  the  river, 
the  highways,  and  the  forestry  of  the  ground  terrain  are 
far  more  recognizable  in  this  deblurred  image.  The 
color  reproduction  of  the  deblurred  image  is  spectacu¬ 
larly  faithful,  and  coherent  artifact  noise  is  virtually 
nonexistent.  There  is,  however,  some  degree  of  color 
deviation  inherently  existing  in  the  deblurred  image. 
These  are  primarily  due  to  chromatic  aberration  and  the 
antireflectance  coating  of  the  transform  lenses.  Nev¬ 
ertheless,  these  two  drawbacks  can  be  overcome  by 
utilizing  good-quality  achromatic  transform  lenses.  A 
research  program  is  currently  underway  to  investigate 
this  effect.  Further  improvement  of  the  deblurring  can 
also  lie  accomplished  by  utilizing  a  blazed  grating  for 
high  diffraction  efficiency  and  a  broader  spatial  band¬ 
width  of  the  deblurred  filter  for  a  higher  degree  of  de¬ 
murring.  These  two  problems  are  also  under  current 
research. 


phase  filter  is  synthesized  by  optical  coating  techniques, 
while  the  fan-shaped  amplitude  filler  is  obtained  by  a 
1-1)  coherent  pnkessing  technique. 

By  comparison  of  the  results  obtained  by  the  broad¬ 
band  image  dehlurring  with  the  narrow  spectral  hand 
and  coherent  techniques,  we  have  seen  that  the  results 
obtained  by  the  broadband  deblurring  offer  a  higher 
image  quality.  We  have  also  shown  fhat  the  broadband 
dehlurring  technique  is  very  suitable  for  color  image 
deblurring.  We  have  provided  several  color  image  de- 
blurring  results  obtained  by  the  broadband  dehlurring 
technique.  From  these  color  deblurred  images  we  have 
seen  that  the  fidelity  of  the  color  reproduct  ion  is  very 
high  ami  the  quality  of  deblurred  image  is  rather  good. 
Although  there  is  some  degree  of  color  blur  due  to 
chromat  ic  aberration  of  the  transform  lenses,  it  can  be 
eliminated  b$-  utilizing  higher-quality  achromatic 
transform  lenses. 

Further  improvements  of  the  dehlurring  can  also  be 
obtained  by  utilizing  a  blazed  grating  to  achieve  a  higher 
smeared  spectral  diffraction  efficiency  so  that  a  wider 
spatial  band  dehlurring  filter  can  he  used  to  achieve  a 
higher  degree  of  dehlurring.  Finally,  we  would  like  to 
point  out  that  the  utilization  of  higher-quality  achro¬ 
matic  transform  lenses  and  a  blazed  grating  for  the 
broadband  image  dehlurring  technique  is  currently 
under  investigation. 


We  acknowledge  the  support  of  the  U.S.  Air  Force 
Office  of  Scientific  Research  grant  AFOSR-81-0148. 


VI.  Summary 

We  have  shown  a  broadband  color  image  dehlurring 
technique  utilizing  a  white-light  source.  'I'his  broad 
spatial  band  dehlurring  technique  utilized  a  grat  ing  base 
method  to  obtain  a  dispersed  smeared  image  spectra  in 
the  Fourier  plane  so  that  the  deblurring  operation  can 
he  taken  placed  in  complex  amplitude  for  the  entire 
visible  wavelengths.  To  perform  this  complex  ampli¬ 
tude  dehlurring  for  the  entire  spectral  band  of  t  he  light 
source,  we  have  shown  that  a  Ian-type  dehlurring  filter 
to  compensate  the  scale  variat  ion  of  the  smeared  signal 
spectra  due  to  wavelength  dispersion  can  be  utilized. 
To  alleviate  the  low  transmission  efficiency  of  (lie  de¬ 
blurring  filter,  we  synthesized  the  dehlurring  filter  with 
the  combination  of  a  broadband  phase  filter  and  a 
fan-shaped  amplitude  filter.  The  broad  spectral  band 
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WHITE-LIGHT  DENSITY  PSEUDOCOLOR  ENCODING 

WITH  THREE  PRIMARY  COLORS 

F.  T.  S.  YU.  X  X  CHF.N.  T.  H  CHAO 
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Traitement  des  images  Optical  processing 

Pseudocoukurt  Pseudocolors 


SUMMARY  :  A  technique  of  generating  density  pseudocolor 
encoding  with  three  primary  colors  using  a  white-light  optical 
processor  is  described.  Spatial  encodings  are  made  with  positive, 
negative,  and  product  of  positive  and  negative  photographic- 
image  transparencies,  and  the  pseudocoloring  is  obtained  by  color 
filtering  of  the  smeared  Fourier  spectra.  The  technique  is  simple, 
versatile,  and  economical  to  operate  that  may  offer  a  wide  range 
of  applications.  Since  coherent  sources  are  not  utilized,  the  color 
coded  image  is  free  from  coherent  artifact  noise.  Experimental 
demonstrations  of  this  density  pseudocolor  encoder  are  provided. 


Images  codecs  par  trois  couleurs  primaires 

RtSUMfi  :  On  decrit  une  nkthode  de  codagc  pseudocolore  avec 
trois  couleurs  primaires  en  ulilisant  un  proeexseur  opliquc  en 
lumicre  blanche.  Le  codagc  pscudocolore  esl  fait  en  fill  rant  par  des 
icrans  colores  les  spectres  de  transparents  pholographiques. 

Le  procede  est  simple  el  il  permet  des  applications  varices  dans  dc 
nombreux  domaines.  L'eclairagc  ilanl  incoherent,  on  ivite  les 
difauts  bien  connus  des  images  en  lumicre  cohercntc 
Quelques  exemples  montrent  Tellicacile  de  la  methode 


Most  of  the  optical  images  obtained  in  various 
scientific  applications  and  usually  gray-level  density 
images.  For  example,  scanning  ejectron  microscopic 
images,  multispcclral  band  aerial  photographic  ima¬ 
ges.  x-ray  transparencies,  etc.  However,  humans  can 
perceive  in  color  better  than  gray-level  variations.  In 
other  words,  a  color  coded  image  can  provide  a 
greater  abilily  in  visual  discrimination. 

In  current  practice,  most  of  the  pscttdocolorings 
are  performed  by  digital  computer  technique  [1].  If 
the  images  arc  initially  digitized,  the  computer  tech¬ 
nique  may  be  a  logical  choice.  However,  for  conti¬ 
nuous  tone  images,  optical  color  encoding  tech¬ 
nique  12]  would  be  more  advantageous  for  at  least 
three  major  reasons  :  first,  the  technique  in  principle 
can  preserve  the  spatial  frequency  resolution  of  the 
image  to  be  color  coded;  second,  the  optical  system 
is  generally  easy  and  economical  to  operate;  third, 
the  cost  of  an  optical  pseudocolor  encoder  is  gene¬ 
rally  less  expensive  as  compared  with  the  digital 
counterpart. 

Density  pseudocolor  encoding  by  halftone  screen 
implementation  with  a  coherent  optical  processor 
was  first  reported  by  Liu  and  Goodman  [3],  and  later 
with  a  white-light  processor  by  Tai,  Yu  and  Chen  [4], 
Although  good  results  have  been  subsequently  report¬ 


ed,  however  there  is  a  spatial  resolution  loss  with 
the  half-tone  technique  and  number  of  discrete  lines 
due  to  sampling  are  generally  present  in  the  color- 
coded  image.  A  technique  of  density  pseudocoloring 
through  contrast  reversal  was  recently  reported  by 
Santamaria  et  al.  [5].  Although  this  technique  offers 
the  advantage  over  the  halftone  technique,  the  optical 
system  is  more  elaborate  and  il  requires  both  inco¬ 
herent  and  coherence  sources.  Since  the  coherence 
source  is  utilized,  the  coherence  artifact  noise  is 
unavoidable. 

More  recently  a  density  pseudocoloi  encoding 
using  a  while-light  processing  technique  was  reported 
by  Chao,  Zhuang  and  Yu  ]6).  This  technique  offers 
the  advantages  of  coherence  noise  reduction,  no 
apparent  resolution  loss,  versatility  and  simplicity  of 
system  operation,  and  low  cost  of  pseudocoloring 
Although  excellent  results  have  been  reported  by  this 
technique,  however  pseudocolor  encoding  is  primarily 
obtained  by  means  two  primary  colors.  We  shall,  in 
this  paper,  extend  this  white-light  pseudocolor  encod¬ 
ing  technique  for  three  primary  colors.  There  is 
however  one  disadvantage  of  this  proposed  technique, 
it  is  still  not  a  real-time  pseudocolor  encoding  tech¬ 
nique. 

We  shall  now  describe  a  simple  white-light  density 


F.  T.  S.  Yu.  X.  X.  Chi  n.  T.  II.  Chau 


J.  Optics  (Paris),  1984.  vol.  15,  n"  2 


pseudocolor  encoding  technique  for  three  primary 
colors.  We  assume  that  a  x-ray  transparency  (called 
a  positive  image)  is  provided  for  pscudocoloring.  By 
contact  printing  process,  a  negative  x-ray  image 
transparency  is  made.  Let  us  now  describe  a  spatial 
encoding  technique  to  obtain  a  threc-gray-lcvcl- 
image  encoding  transparency  for  the  pscudocoloring. 
The  spatial  encoding  is  performed  by  respectively 
sampling  the  positive,  negative,  and  the  combination 
of  both  (i  e„  the  product  of  positive  and  negative) 
transparencies  onto  a  black-and-white  photographic 
film,  with  specific  sampling  grating  frequencies  orient¬ 
ed  at  specific  azimuthal  directions.  To  avoid  the 
Moire  fringe  pattern,  we  shall  sample  these  three 
images  in  orthogonal  direction  with  different  specific 
sampling  frequencies,  as  proposed  in  figure  I.  As  an 
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illustration,  the  intensity  transmittance  of  the  encoded 
film  can  be  written  as 

(I)  7  (x,  y)  =  K  {  7  ,(x.  y)[l  +  sgn  (cos /»,  >)] 

+  7  2(x,  y)  [1  +  sgn  (cosp2x)J 
+  T 3(x.  >) (1  +  sgn  (cosp3  x)]  }'y, 

where  K  is  an  appropriate  proportionality  constant, 
T,.  T j,  and  T,  are  the  positive,  the  negative,  and 
the  product  image  exposures,  />,.  p2,  and  p3  are  the 
respective  earner  spatial  frequencies,  (x.  y)  is  the 
spatial  frequency  coordinate  system  of  the  encoded 
film.  >■  is  the  film  gamma,  and 


sgn  (cos  v )  £ 


1,  cos  x  >  0 , 
-  1,  cos  x  <  0  . 


We  shall  now  bleach  the  encoded  transparency  to 
obtain  a  surface  relief  phase  object  [7,  X).  We  assume 
that  the  bleached  transparency  is  encoded  in  the 
linear  region  of  the  diffraction  efficiency  D  versus 
log-exposure  E  curve  (8).  Thus,  the  amplitude  trans¬ 
mittance  of  the  bleached  transparency  can  be  written 
as 

(3)  r(x,  y)  =  exp[i$(x,  y)]  , 

where  (f>(x.  y)  represents  the  phase  delay  distribution, 
which  is  proportional  to  the  exposure  of  the  encoded 
film  [9],  such  as 

(4)  y)  =  M  {  7  ,(x.  y)  l*  +  sgn  (cos/>,  y)J 

+  7  ,(.y.  y)  ( I  +  sgn  (cos/>2  x)] 

+  7  jf.v,  y)  11  +  sgn  (cos  ps  x)J  } , 


E»i«nded  White-Light 
Sourct 


Negolivi  J  ^-Posii*** 
TrohipOfetKt  /  Trah*00f«nc» 


®  Blue  Co'or  F iit*» 


Groen  Color 
^  F'l*«r 


/y* 
*  ; 

LI 

■f  J-  . 

S'  Output 

Pi 

1 

\  N.  N  j  *\ 
- 

1 

f 

1  Ip, 

l;Ki.  2.  A  whitv-lighf  p.wuthn  otor  vntothr 


where  M  is  an  appropriate  proportionality  constant 
If  we  place  this  bleached  encoded  film  at  the  input 
plane  Px  of  a  white-light  optical  processor,  as  illus¬ 
trated  in  figure  2  then  the  complex  light  distribution 
due  to  f(x,  y),  for  every  A,  at  the  spatial  frequency 
plane  P2  can  be  determined  by  the  following  Fourier 
transformation  : 


(5)  S(ar.  P :  A)  =  J  j  f(x,  y)  expj^  -  i  jj  (ax  +  0y)J  dx  dy  =  j  j  exp[id»(x.  v )]  expj^  -  i  (ax  +  fly)  j  dx  dy 
By  expanding  f(x,  y)  into  an  exponential  series,  Eq  (5)  can  be  written  as 

(6)  S(a.  p ;  A)  =  j  j  j  )  +  y)  +  j  [i< M*.  y)]2  +  -  |  expj^  -  i  (ax  +  fly  )j  dx  dy  • 

By  substituting  Eq.  (4)  into  Eq.  (6)  and  retaining  the  first-order  terms  and  the  first-order  convolution  terms,  we 
have  : 
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where  T,.  T 2  and  T 3  arc  the  Fourier  transforms  of  T,.  T"2  and  T 3  respectively.  *  denotes  the  convolution  opera¬ 
tion.  and  the  proportional  constants  have  been  neglected  for  simplicity.  We  note  that,  the  last  cross  product  term 
of  Iq.  (7)  would  introduce  a  Moire  fringe  pattern,  which  in  the  same  sampling  direction  of  />2  and  />,.  Nevertheless, 
all  of  these  cross  product  terms  can  he  properly  masked  out  at  the  Fourier  plane.  Thus  hy  proper  color  filtering 
the  first-order  smeared  Fourier  spectra,  as  shown  in/igwiv  .f,  a  Moire  free  pseudocolor  coded  image  can  be  obtained 
at  the  output  plane  P3.  The  corresponding  complex  light  field  immediately  behind  the  Fourier  plane  would  be 


(8) 
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where  Xr.  Xh  and  arc  the  respective  red,  blue,  and 
green  color  wavelengths.  At  the  output  image  plane, 
the  pseudocolor  coded  image  irradiance  is  therefore, 

(9)  / (x,  y)  =  T j,(x.  y)  +  T |„(x,  y)  +  T f,(x,  y) , 

which  is  a  superposition  of  three  primary  color 
encoded  images,  where  7  ,r.  Tlb.  and  Tit  are  the 
red,  blue,  and  green  amplitude  distributions  of  the 
three  spatially  encoded  images.  Thus  a  Moire  free 
color  coded  image  can  be  obtained  at  the  output 
plane. 

In  our  experiment,  we  utilized  two  sinusoidal 
sampling  gratings,  one  is  26.7  lines/mm  and  the 
other  is  40  lines/mm  for  the  spatial  encodings.  The 
encoding  transparency  was  made  by  Kodak  tech¬ 
nical  pan  film  2415.  The  advantage  of  using  Kodak 


2415  film  is  that  it  is  a  low  contrast  film  with  a  rela¬ 
tively  fiat  spectral  response.  The  plot  of  diffraction 
efficiency  versus  log-exposure  for  Kodak  2415  film 
at  40  lincs/mm  sampling  frequency  is  shown  in 
figure  4.  From  this  figure,  we  see  that  the  bleached 
encoded  films  offer  a  higher  diffraction  efficiency,  the 
optimum  value  occurs  at  exposures  8.50  x  I0~*  mcs. 
With  reference  to  this  optimum  exposure,  it  is  possible 
to  optimize  the  encoding  process  in  the  following  : 
first,  preexposuring  the  film  beyond  the  shoulder 
region ;  second,  subdividing  the  remaining  exposure 
into  three  regions  by  taking  the  account  of  the  trans- 
mitlcnt  exposures  of  (he  three  encoded  images. 


Log  Exposure! mes) 


Ik;.  4  Dill  rod  ion  i  ffu  n/ny  versus  log  c\ junior  /»/<</  of  a  iypitol 

spat  tally •  cm  oih  '(/film 


In  pseudocoloring,  we  utilize  the  Kodak  primary 
color  filters  of  No.  25,  47B  and  58  in  the  Fourier 
plane,  as  shown  in  figure  3.  A  xenon-arc  lamp  is  used 
as  extended  white-light  source  for  pseud  lor  encod¬ 
ed  of  figure  2. 

For  experimental  demonstrations,  we  would  first 
provide  a  gray-level  x-ray  picture  of  a  front  view 
female  pelvis,  as  shown  in  figure  5(a).  Figure  5(b) 
shows  the  color  encoded  image  obtained  by  the 
white-light  pseudocolor  encoder,  in  (his  color  coded 
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Klli  5.  —  (a)  A  gruv-lcvcl  v-niy  picture  (b).  A  density  pseudocolor 
cmlnl  pit  lure  Positive  image  it  coded  m  ml.  negatin'  image  is  axled 
in  blue.  mid  the  prodtn  I  iinagt  r,  axial  in  gra  il  It)  A  reversal  color 
axial  initigc  n/  figure  5  (b)  Puvilne  image  is  axial  in  blue,  negative 
image  is  t  oded  m  ral.  and  the  prnilui  I  image  i.\  again  coded  in  green 
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image,  the  positive  image  is  encoded  in  red,  the 
negative  image  is  encoded  in  blue,  and  the  product 
transmittance  of  the  positive  and  negative  images  is 
encoded  in  green.  From  this  figure,  we  see  that  a 
broad  range  of  pseudocolor  encoded  density  can  be 
perceived,  and  the  color-coded  image  appears  free 
from  coherent  artifact  noise  and  Moire  fringes.  The 
patient  of  the  x-ray  picture  appears  to  have  suffered 
from  a  surgical  operation,  where  the  bore  and  point 
hardwares  can  be  seen.  A  section  of  the  bone,  between 
the  sacroiliac  joint  and  spinal  column,  has  been 
removed. 

It  may  also  be  interesting  to  note  that  a  reversal 
of  the  color  encoding  can  be  easily  obtained  by  this 
whitc-ligh'  pseudocolor  encoding  technique,  as  shown 
in  figure  5(c).  In  this  figure,  the  positive. image  is 
encoded  in  blue,  the  negative  image  is  encoded  in 
red,  and  the  product  transmittance  is  encoded  in 
green.  From  this  figure,  again  we  see  that  a  broad 
range  density  pseudocolor  encoded  image  with  diffe¬ 
rent  color  textures  can  be  obtained.  For  example  the 
air  pockets  in  the  colon  of  the  patient  in  figure  Sic) 
is  much  easier  to  be  identified.  In  concluding  this 
paper,  we  note  that  a  wide  variety  of  pseudocolor 
encoded  images  can  easily  be  obtained  by  simply 
alternating  the  color  filters  in  the  Fourier  plane  of 
the  while-light  processor.  The  proposed  white-light 
pseudocolor  encoder  is  economical  and  easy  to 
operate,  which  may  offer  many  practical  applications. 
There  is.  however,  one  disadvantage  of  this  technique, 
it  requires  a  spatial  encoding  process.  Therefore, 
it  is  not  a  real-time  pseudocolor  encoding  method. 

We  acknowledge  the  support  of  the  U.S.  Air  Force 
Office  of  Scientific  Research  Grant  AFOSR-8I-OI48. 
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Coherence  measurement  of  a  grating-based  white-light 
optical  signal  processor 


F.  T.  S.  Yu,  F.  K.  Hsu,  and  T.  H,  Chao 


A  two-licnm  inlcri'rronrr  technique  lor  coherence  iiienxiirenienl  in  I  In*  Fourier  plane  of  n  grating- liaxetl 
white-light  optical  signal  processor  is  presented.  The  visibility  measurement  is  ohlnined  with  n  scanning 
pholomeler  Hi  (lie  output  plane  of  the  processor.  The  degree  of  coherence  ns  n  function  of  slit  separation, 
due  to  source  sue,  input  object  size,  and  the  spatial  frequency  of  the  sampling  gral  ing.  is  plolted.  The  re¬ 
sults  show  that  the  degree  of  coherence  increases  as  the  spatial  frequency  of  the  sampling  grating  increases. 
However,  this  improvement  of  coherence  is  somewhat  more  effective  in  the  direction  perpendicular  to  the 
light  dispersion.  Thus,  this  white-light  signal  processing  technique  is  more  effective  in  one  dimension. 
Since  this  white-light  processor  is  capable  of  processing  the  signal  with  the  entire  spectral  band  of  the  light 
source,  it  is  very  suitable  for  color  image  processing. 


Introduction 


In  the  19.30s,  work  by  Van  Ciltert1  and  later  by  Zer- 
nike-  drew  attention  to  the  study  of  partial  coherence. 
They  have  shown  that  the  spatial  coherence  and  the 
intensity  distribution  of  the  light  source  form  a  Fourier 
transform  relationship.  More  recently,  the  work  of 
Thompson*' s  has  demonstrated  a  two-beam  interfer¬ 
ence  technique  to  measure  the  degree  of  partial  coher¬ 
ence.  He  has  shown  that,  under  quasi-momiehromalic 
illumination,  the  degree  of  spatial  coherence  is  depen¬ 
dent  on  the  source  si/e  and  the  distance  between  two 
arbitrary  points.  The  degree  of  temporal  coherence  is 
however  dependent  on  the  spectral  bandwidth  of  the 
light  source.  He  has  also  illust  rated  several  coherence 
measurements  that  are  consistent  with  the  Van  Cittert 
mode  in  the  Fourier  predictions. 

We  recently  proposed  a  grating-based  white-light 
optical  signal  processor.*’  H  We  have  shown  that  the 
white-light  processor  is  capable  of  processing  the  in¬ 
formation  in  complex  amplitude  like  a  coherent  pro¬ 
cessor,  and  at  the  same  time  it  suppresses  the  coherent 
artifact  noise  like  an  incoherent  processor.  In  other 
words,  this  white-light  processor  is  operating  in  a  par¬ 


tially  coherent  mode  in  the  Fourier  plane  as  a  partially 
coherent  processor  instead  of  an  incoherent  pro¬ 
cessor. 

In  this  paper  we  shall  describe  a  dual-beam  technique 
for  coherence  measurement  in  the  Fourier  plane.  We 
shall  show  that  a  high  degree  of  coherence  can  be  ob¬ 
tained  in  the  spatial  frequency  plane  such  that  the  sig¬ 
nal  can  be  processed  in  complex  amplitude  for  the  entire 
spectral  hand  of  t  he  white-light  source. 


II.  White-Light  Optical  Processing  Technique 

With  reference  to  the  schematic  diagram  of  Fig.  1,  the 
white-light  optical  processor  is  similar  to  that  of  a  co¬ 
herent  opt  ical  processing  system,  except  lor  the  use  of 
an  extended  white-light  source,  source  encoding  mask, 
signal  sampling  grating,  and  achromatic  transform 
lenses.  By  the  Wolf  partial  coherence  theory,9  the 
output  intensity  distribution  can  he  obtained  by  the 
following  integrating  equations1": 


Hx’y)  *  jy  7(xn.yo) 
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where  tUoO'o)  *s  the  intensity  distribution  of  the  source 
encoder  at  the  source  plane  UovVo).  S(a,f})  is  the  Fourier 
spectrum  of  the  input  signal  s(x,y),  Pn  is  the  angular 
spatial  frequency  of  the  grating,  /  is  the  focal  length  of 
the  achromatic  transform  lenses,  // (o,d)  is  the  complex 
spatial  filter  in  the  Fourier  plane,  \/,\h  are  the  lower 
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Fig.  I.  Grating-based  white-light  opti<~il  signal  processor:  -yCrcyVn), 
source  encoding  mask;  L,  achromatic  transform  lenses. 


and  upper  wavelength  limits  of  the  light  sotiree,  nod 
(u,/f)  is  the  spatial  coordinate  system  of  the  Fourier 
plane. 

In  optical  signal  processing,  a  set  of  N  narrow  discrete 
spectral  band  filters  H„  (<v„  ,0„ )  over  the  rainbow  color 
of  Fourier  spectra  are  used  in  the  Fourier  plane.  Each 
of  the  spectral  band  filters  is  limited  by  a  finite  spectral 
bandwidth  that  can  be  approximated  by  the  following 
equation7-8: 

4A  p 

AX  as - X„ ,  po  »  A p.  (2) 

Po 

where  A  p  is  the  angular  spatial  frequency  limit  of  the 
input  signal  in  the  x  direction,  p0  is  the  angular  spatial 
frequency  of  the  sampling  grating,  and  A„  is  the  center 
wavelength  of  the  filter  Hn  («n  ,0n ).  Since  the  narrow 
spectral  band  filtered  signals  are  mutually  incoherent, 
the  output  intensity  distribution  of  Eq.  (1)  can  be 
written  as 


/U'.Jt')=*  I  )  ==  V  j  If  •>(*„, v„) 

»i“l  m»1  AW?  »/»/-•» 

•  j  jj  n(v„  +  «  -  ~  PikVo  +  |t| II „  (<>„  ,rf„ ) 
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With  reference  to  Eq.  (1)  or  (3),  we  see  that  the 
white-light  signal  processor  is  indeed  capable  of  pro¬ 
cessing  the  signal  in  complex  amplitude  for  the  entire 
spectral  band  of  the  light  source,  and  it  is  very  suit  able 
for  color  image  processing.  From  these  two  equations, 
we  also  see  that  the  degree  of  coherence  is  governed  by 
the  source  encoding  mask  7(x(l,yo),  the  spatial  fre¬ 
quency  of  the  sampling  grating  p(>.  and  the  spatial  fre¬ 
quency  limit  of  I  lie  input  signal  Ap.  In  oilier  words,  the 
source  encoding  function  is  to  improve  the  degree  of 
sputial  coherence  ul  the  input  plunc  and  the  signal 
sampling  grating  is  to  increase  the  degree  of  temporal 
coherence  in  the  Fourier  plane,  so  that  the  signal  pro¬ 
cessing  can  be  carried  out  in  complex  amplitude  at  each 
narrow  spectral  band  filler. 

It  is  not  hard  to  see  from  Fig.  1  that  the  signal  spec¬ 
trum  is  dispersed  into  rainbow  color  along  the  a  axis  in 
the  Fourier  plane.  The  improvement  of  the  degree  of 
coherence  is  expected  to  be  more  effective  in  the  0  di¬ 
rection  than  in  the  cv  direction.  In  devising  a  coherence 
measurement  scheme  for  this  white-light  optical  pro- 
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lo>  Input  Plone  ipi  Fourier  Plone 

Fig.  2.  Visibility  measurement  along  the  fi directions:  l),  input  slit 
width;  d,  mean  slit  serration;  Xr>.  mean  wavelength  of  the  light  source; 
Pu.  angular  spatial  fre«|iicnov  n  f  I  lie  iiniii|iling  grilling. 

cessor,  we  propose  n  dunl-liciun  interference  technique 
for  coherence  measurement  in  the  Fourier  plane.  For 
simplicity,  we  would  however  measure  the  degree  of 
coherence  in  the  0  and  in  the  a  axes  independently. 

III.  Coherence  Measurement  Technique 

We  shaH  now  describe  a  visibility  measurement 
technique  to  determine  the  degree  of  coherence  in  the 
0  direction  in  the  Fourier  plane.  For  simplicity,  we 
shall  utilize  a  narrow  slit  as  a  1-D  object  at  the  input 
plane  as  shown  in  Fig.  2(a).  The  complex  light  distri¬ 
bution  in  the  Fourier  piane  would  be 

£(nr,0;X)  «  c|sinc  /sjj  •  6  |or  -  ^  •  (4) 

where  C  represents  an  appropriate  complex  constant, 
D  is  the  slit  width,  f  is  the  focal  length  of  the  achromatic 
transform  lens,  p0  is  the  angular  spat  ini  frequency  of  the 
phase  sampling  grating,  and  *  denotes  the  convolution 
operation. 

It  is  clear  that  Eq.  4  describes  a  fan-shaped  smeared 
Fourier  spectrum  of  the  input  slit,  for  which  the  scale 
of  the  sine  fnctor  increases  ns  n  function  of  the  wave¬ 
length  \  of  the  light  source  and  decreases  as  the  size  of 
the  object  D  (i.e.,  slit  width)  .ncreases.  To  increase  the 
efficiency  of  the  coherence  measurement  along  the  0 
direction,  we  would  use  n  pair  of  slanted  narrow  slits  nt 
the  smeared  Fourier  spectra,  as  illustrated  in  Fig.  2(b). 
The  angle  of  inclination  of  this  pair  of  slits  should  he 
ud justed  with  the  separation  of  the  slits,  such  as, 

0  ■  inn  ■■  .  tr») 

An//>t» 

where  (l  is  the  menn  sepimilion  of  l ho  slits*  A,»  is  the 
mean  wavelength  of  the  light  source,  and  p()  is  the  an¬ 
gular  spatial  frequency  of  the  sampling  grating.  The 
filtering  function  of  this  pair  of  slanted  slits  can  he 
written  as 

+  ,6, 

The  output  intensity  distribution  can  be  shown  as 

/(*',>’')“  K[1  +  cos(2Tdx')),  (7) 

where  K  is  an  appropriate  proportionality  constant. 
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Kquat  ion  (7)  represents  nn  achromatic  fringe  pul  torn 
due  lo  this  entire  spectral  hand  of  the  light  source. 

Strictly  speaking,  all  practical  white- light  sources  are 
extended  sources.  For  simplicity  of  illustration,  we 
assume  that  an  extended  square  source  is  used  in  this 
white-light  optical  processor.  Thus  the  output  inten¬ 
sity  distribution  can  be  written  as 

/(tv>  -  xr  JX’ rec*  H'H?) 

’  I JX-  sinc  [ly  ^  +  >n)]  *  4(°  +  *°  -  ^  Po.fl) 

•W(or,0)exp|-i^y(jr'or  +  .v'0)  dodtfj  dxodyodX,  (8) 
where  H(a,0)  is  given  by  Eq.  (6), 

rod  (— ]  -  (*'  jXnj  “  a/^‘ 

In/  lo.  |.v„|  >  n/2, 

and  a  is  the  dimension  of  the  extended  light  source. 
From  the  above  equation  we  see  that  the  visibility  (i.e., 
degree  of  coherence)  is  dependent  on  the  source  si/.eri, 
the  object  size  I)  (i.e.,  slit  width),  and  the  angular  spatial 
frequency  p(t  of  the  sampling  grating. 
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To  investigate  the  degree  of  coherence  variation  in 
the  nr  direction,  again  we  insert  a  narrow  slit  aperture 
ns  nn  input  object  but  parallel  to  the  sampling  direction 
of  the  sampling  grating,  as  shown  in  Fig.  3(a).  The 
smeared  Fourier  spectra  can  be  shown  as 

£(«,/);  A I  -  <«j  ♦  *|«r  po.dj  •  (9) 

which  descriltcs  n  narrow  smeared  rninlmw  color  s|>ectrn 
along  the  o  axis.  For  the  visibility  measurement,  a  pair 
of  narrow  slits  is  inserted  at  the  Fourier  plane  perpen¬ 
dicular  to  the  «  axis  and  centered  at «  *  Aa/p0/( 2x),  as 
shown  in  Fig.  3(b),  where  A0  is  the  center  wavelength  of 
t  he  light  source  and  { is  the  focal  lengt  h.  The  filtering 
function  of  this  pair  of  slits  cun  he  written  us 

H(n,0)  -  -  -  -  ^p0,^J  +  a|,v  +  £-  T^Pu.dj  •  (10) 

where  li  is  the  separation  between  the  slits.  Again  with 
the  assumption  of  a  square  light  source,  the  output  in¬ 
tensity  distribution  can  be  shown  as 

•  |  XT- Binc  [^7 {a + io)J  *  4*  -  77 p"-  * + 

12  IT  1  2 

-i  —  (x'or  +  y'/SlJdmffl  dxodvodA,  (11) 

where  H («,/?)  is  defined  in  Eq.  (10).  From  this  equa¬ 
tion,  again  we  see  that  the  degree  of  coherence  in  the  « 
direction  is  dependent  on  the  source  size  a,  the  object 
size  D,  and  the  angular  spatial  frequency  p0  of  the 
sampling  grating. 

IV.  Experimental  Results 

The  optical  setup  for  the  measurement  of  the  degree 
of  coherence  in  the  Fourier  plane  or  a  grating-bused 
white-light  option)  processor  is  shown  in  Fig.  <1.  'I'his 
setup  utilizes  the  principle  of  a  dual-beam  interference 


Fig.  4.  Optical  setup  for  the  coherence  measure¬ 
ment  along  the  0  direction:  ■y(.to.yn),  source  en¬ 
coding  mask;  a,  source  size;  I).  input  slit  size;  L, 
achromatic  transform  lenses;  H(n,0)  pair  of  slant 
slits;  l‘M,  photometer;  OSC,  oscilloscope. 
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Fig.  5.  1’lnlsof  Hegrer  of  rohrrpncr  nlong  I  ho  fi  Hirorl  ion  ns  n  func¬ 
tion  nf  nii'nit  slit  srpnrnl ion  1/  lor  various  vnlursol  source  si/cti.  f,i 
*  d/(  W>.  the  mean  slil  separation  in  spatial  frequency. 


technique  for  coherence  measurement.  The  output 
interference  fringe  pattern  can  be  traced  by  a  linear 
scanning  photometer  and  displayed  on  an  oscilloscope 
for  the  visibility  measurement.  In  the  experiment,  the 
photometer  is  made  by  mounting  a  photomultiplier  on 
the  top  of  a  motor-driven  linear  translator.  Since  it  is 
a  1  -D  fringe  pattern,  a  narrow  slit  can  be  utilized  at  the 
input  end  of  the  photomultiplier  for  the  visibility 
measurement,  i.e., 


where  /„„„  and  are  the  maximum  and  minimum 
intensities  of  I  In*  hinges.  Needless  In  soy  that  the  vis¬ 
ibility  of  the  fringes  is  in  fact  a  measure  of  t  he  degree  of 
coherence  measurement,'*0  i.e.,  V  =  | > | ,  where  7  is  the 
complex  degree  of  coherence. 

In  the  following  we  shall  illustrate  (he  visibility  (i.e., 
coherence)  measurement  in  the  ft  and  «  directions  in  the 
Fourier  plane. 

A.  Case  I:  Coherence  Measurement  in  the  ft 
Direction 

We  shall  now  describe  the  coherence  measurement 
in  the  ft  direction  as  illustrated  in  Fig.  4.  We  shall  show 
that  the  visibility  varies  as  functions  of  mean  separation 
d  of  the  pair  of  slanted  slits,  source  size  a,  object  size  D 
(i.e.,  slit  width),  and  spatial  frequencies  of  the  sampling 
grating.  Figure  5  shows  the  variation  of  the  degree  of 
coherence  as  a  funct  ion  of  main  separat  ion  d  for  various 
values  of  source  sizes  a.  From  this  figure  we  sec  that 
the  degree  of  coherence  decreases  as  the  separation  d 
increases.  Further  increase  in  d  increases  the  reap¬ 
pearance  of  the  coherence.  Still  further  increase  in  d 
causes  the  repeated  fluctuation  of  visibility.  In  this 
figure,  we  also  see  that  the  degree  of  coherence  increases 
as  the  source  size  a  decreases.  There  is  an  interesting 
phenomenon  in  this  coherence  measurement.  We  see 
that  as  the  source  size  a  decreases  further,  the  reap- 


|N'iirmHi’  of  the  visibility  is  higher.  I  Inwovrr  tin*  ovrrnll 
intensity  of  the  smeared  Fourier  spectra  decreases. 
This  phenomenon  is  primarily  due  to  the  finite  object 
size  under  a  uniform  source  size  illumination.  Fur¬ 
thermore,  if  the  source  size  a  is  further  increased,  for 
example,  exceeding  1.0  mm  in  Fig.  5,  the  decrease  in 
.  oherence  due  to  the  source  size  is  not  apparent.  This 
is  primarily  due  to  a  comparable  broader  object  size  D 
(e.g.,  D  =  0.6  mm)  compared  to  a  narrower  sine  factor 
derived  from  source  size  a.  Nevertheless,  if  the  input 
slit  size  further  decreases,  the  changes  in  degree  of  co¬ 
herence  for  large  source  sizes  can  he  seen. 

We  shall  now  provide  a  set  of  output  fringe  patterns 
with  a  set  of  normalized  scanned  photometer  traces  as 
shown  in  Fig.  6.  The  fringe  pul  terns  of  Figs.  6(ji)  (<•) 
were  taken  at  visibilities  of  0.88,  0.68,  ami  0.08  that 
corres|x>nd  to  the  mean  separations  d  -  0.116, 0.72,  and 
1.41  mm  at  (minis 0,  />,  and  c  indicated  on  the  main  Inin* 
of  the  plot  a  -  0.4  mm  in  Fig.  5.  Figures  6(d)  and  (e) 
were  taken  on  the  second  lobe  of  the  visibility  reap¬ 
pearance  that  correspond  to  points  d  and  e  in  Fig.  5. 
The  degrees  of  coherence  at  these  two  points  are  0.42 
and  0.15,  respectively.  And  the  corresponding  mean  slit 
separations  are  2.16  and  2.52  mm.  Figure  6(f)  was  taken 
on  the  third  lobe  visibility  at  point  f.  The  degree  of 
coherence  is  0.3  and  the  mean  slit  separation  is  2.88 
mm. 

Let  us  now  investigate  the  degree  of  coherence  as  a 
function  of  mean  separation  d  for  various  input  object 
sizes  U  (i.e.,  slit  width)  as  plotted  in  Fig.  7.  Again,  we 
see  that  the  degree  of  coherence  decreases  as  d  increases. 
Further  increase  in  d  also  causes  side  lobes  to  reappear. 
In  this  figure  we  also  see  that  the  degree  of  coherence 
increases  as  object  size  D  decreases.  Finally,  Fig.  8 
shows  the  visibility  measure  as  a  function  of  mean 
separation  d  for  two  values  of  sampling  grating 
frequencies.  From  this  figure  we  see  that  the  degree  of 
coherence  is  dramat  ically  improved  in  the  Fourier  plane 
by  the  insertion  of  the  sampling  grating. 

B.  Case  II:  Coherence  Measurement  in  the  « 
Direction 

At  this  stage  we  shall  now  measure  the  degree  of  co¬ 
herence  in  the  <v  direction  in  the  Fourier  plane.  The 
measurement  technique  is  essentially  identical  to  that 
of  Fig.  4,  except  the  input  object  is  replaced  by  a  pair  of 
horizontal  slits  as  shown  in  Fig.  3.  In  coherence  mea¬ 
surement,  we  centered  the  pair  of  slits  at  the  center  of 
the  smeared  Fourier  spectra  corresponding  to  X  =  5461 

A. 

Figure  9  shows  plots  of  degree  of  coherence  as  a 
function  of  slil  separation  Ii  l  or  various  values  of  source 
sizes  0.  From  this  figure  we  see  that  the  degree  of  co¬ 
herence  decreases  as  h  increases.  Further  increase  in 
h  again  causes  the  reappearance  of  the  visibility  side 
lobes.  However,  the  degree  of  coherence  is  generally 
not  affected  by  the  variation  of  the  source  size  0 .  Figure 
10  shows  a  set  of  the  visibility  fringe  patterns  that  we 
have  obtained  in  the  output  image  plane.  This  set  of 
pictures  was  taken  at  points  a,  b,  c,  and  d  as  shown  in 
Fig.  9.  The  corresponding  degrees  of  coherence  are 
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rip.  6.  Samples  of  frinpe  visibility  patterns  at  the  output  plane.  1  he  upper  portion  of  (a)  -(f).  shows  the  frinpe  visibility  patterns.  The  lower 
portion  shows  the  corresponding  photometer  traces.  In  these  experiments,  the  fringe  visibility  and  the  spatial  frequency  were  varied  by  changing 
the  mean  separation  distance  d  between  the  two  slanted  slits,  (a)-(c)  were  obtained  at  points  o,  b,  and  c  on  the  first  lobe  of  Fig.  5.  These 
figures  show  the  decrease  in  fringe  visibility  and  corresponding  increase  in  spatial  frequency  as  the  separation  d  increases,  (d)  and  (e)  were 
obtained  at  ooints  d  and  c  on  the  second  lolie  of  Fig.  f>.  These  two  figures  also  show  the  increase  of  spatial  frequency  as  d  further  increases. 

(f)  was  obtained  at  point  f  on  the  third  lobe  of  Fig.  5. 
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Kig.  7.  Plots  of  degree  of  coherence  along  the  (i  direction  as  a  funct¬ 
ion  of  mean  alit  separation  d  for  various  values  of  input  slit  width 

D. 


Kig.  8.  Plots  of  degree  of  coherence  along  the  d  direction  as  a  func¬ 
tion  of  mean  slit  separation  d  for  two  values  of  sampling  grating  fre¬ 


quency  po. 


Fig.  9.  Plots  of  degree  of  coherence  along  the  or  direction  as  a  func¬ 
tion  of  slit  separation  h  for  various  values  of  source  size  o.  f„  "  h! 
(Xof),  the  corresponding  slit  separation  in  spatial  frequency. 


().7f>,  ()..r>0,  O.M,  und  0.05.  The  rcs|>ective  separations 
are  h  =  0.36, 0.72, 1.44,  and  2.16  mm. 

We  shall  now  plot  the  degree  of  coherence  as  a  func- 
tion  of  h  for  various  values  of  object  sizes  I)  (i.c..  slit 
widths  D)  as  shown  in  Fig.  11.  From  these  plots  we  see 
that  the  degree  of  coherence  decreases  as  the  object  size 
I)  increases.  Figure  12  shows  the  variation  of  coherence 
due  to  spatial  frequency  of  the  sampling  grating  as  a 
function  of  slit  separation  h.  From  this  figure  we  see 
that  higher  degree  of  coherence  is  achievable  with  the 
insertion  of  a  high  spatial  frequency  grating  in  the  input 
plane. 

We  shall  now  briefly  discuss  the  overall  effect  of  co¬ 
herence  in  the  ( at,0 )  spatial  frequency  plane.  By  com¬ 
paring  t  he  visibility  measurement  of  Figs.  5  and  0.  we 
see  t  hat  t  he  degree  of  coherence  substant  ially  increases 
in  t  he  0  direct  ion  as  the  source  size  decreases.  There 
is,  however,  no  significant  improvement  in  the  <r  di¬ 
rection  for  smaller  source  sizes.  Although  both  cases 
show  the  increase  in  coherence  for  smaller  object  sizes, 
however,  the  increase  in  coherence  is  higher  in  the  0 
direction  compared  with  the  a  direction,  as  shown  in 
Figs.  7  and  11.  With  reference  to  the  plots  of  Figs.  8  and 
12,  both  cases  show  significant  improvement  in  the 
degree  of  coherence  with  the  insertion  of  a  sampling 
grating.  However  the  improvement  in  coherence  in  the 
0  direction  is  somewhat  higher  than  in  the  a  direction, 
due  primarily  to  overlapping  of  the  smeared  rainbow 
Fourier  spectra. 

To  summarize  these  observations,  we  stress  that  the 
grating-based  white-light  optical  signal  processor  does 
improve  the  degree  of  coherence  in  the  Fourier  plane. 
Although  the  higher  degree  of  coherence  is  obtainable 
in  both  spatial  frequency  directions,  the  coherence 
improvement  in  the  0  direction  is  generally  higher. 
Thus,  the  white-light  optical  signal  processing  tech¬ 
nique  is  generally  more  effective  in  the  0  direction  than 
in  the  or  direction.  We  note  that  this  effect  can  he  seen 
in  a  recent  paper  on  linear  motion  color  image  de- 
blurring.11 


V.  Conclusion 

We  have  devised  a  dunl-benm  interference  technique 
to  measure  the  degree  of  coherence  in  the  Fourier  plnne 
of  a  grating-based  white-light  optical  signal  processor. 
The  effect  of  coherence  variation  due  to  source  size, 
input  object  size,  and  the  spatial  frequency  of  the 
sampling  grating  is  plotted  ns  a  function  of  distance  in 
the  0  and  or  directions  of  the  Fourier  plane.  We  have 
shown  that  the  degree  of  coherence  increases  as  the 
spatial  frequency  of  the  sampling  grating  increases. 
Although  the  improvement  in  degree  of  coherence  in  the 
Fourier  plane  is  quite  evident,  the  improvement  in  the 
0 direction  (i.e.,  the  direction  perpendicular  to  the  light 
dispersion)  is  somewhat  more  effective  than  in  the  « 
direction.  The  results  indicate  that  this  white-light 
optical  signal  processing  technique  is  somewhat  more 
effective  in  the  0  direction  than  in  the  n  direction. 
Nevertheless  the  existence  of  the  high  degree  of  coher¬ 
ence  in  the  Fourier  plane  allows  us  to  process  the  in- 
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Fin.  10.  Samples  of  fringe'  visibility  patterns.  TIm>  upper  |mrtion  of  (a)  (d)  shows  frinji*"  visibility  pall  crus;  (he  lower  imrtion  shows  the  cor¬ 
responding  intensity  profiles,  (a)  and  (h)  were  obtained  at  points  a  and  b  on  the  first  l»l»e  of  Fig. !).  (c)  was  ohlnined  at  point  c  on  the  second 

lolic  of  Fig.  0,  mid  (d)  wns  ohlnined  at  point  d  on  the  third  lolie  of  Fig. !). 


o  J«  'os*  iss«  ITS0  t«  *hA.t(<'"««  t  pjg  j2.  plots  of  degree  of  coherence  along  the  o  direction  as  a 

Fig.  11.  Plots  of  degree  of  coherence  along  the  or  direction  as  a  function  of  slit  separation  h  for  two  values  of  sampling  grating  fre- 
function  of  slit  separation  h  for  various  values  of  object  size  D.  quency  po¬ 
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A  method  of  deconvolving  the  dcfocuscd  color  photographic  images  utilizing  a  white-light  processing  is  described.  In 
the  white-light  processing,  a  diffraction  grating  is  used  to  provide  three  primary  color  light  sources  for  the  color  image 
restoration.  Three  complex  inverse  filters,  for  each  primary  color,  arc  used  in  the  Courier  transform  plane.  Experimental 
demonstrations  of  the  color  image  restoration  of  dcfocuscd  photographic  images  arc  given. 


1 .  Introduction 

Restoration  of  blurred  photographic  images  has 
lung  been  an  interesting  application  in  optical  signal 
processing  [I -8).  We  have  in  previous  papers  (9-1 2J 
presented  a  white-light  processing  technique  for  lin¬ 
early  smeared  image  deblurring.  We  have  shown  that 
the  white -light  image  dcblurring  technique  is  capable 
of  suppressing  the  coherent  artifact  noise  and  is  suit¬ 
able  lor  color  Image  dcblurring.  However  the  results 
that  we  had  obtained  were  primarily  restricted  to 
binned  due  to  lineal  motion. 

In  this  paper,  we  shall  extend  the  white-light  image 
restoration  technique  to  two-dimensional  image  de- 
blurring.  We  shall  utilize  a  whito-light.sourcc  for  the 
restoration  of  the  out-of-fdcuscd  color  photograph 
images.  The  proposed  color  image  restoration  tech¬ 
nique  includes  three  primary  color  sensitive  inverse 
filters,  as  shown  in  fig.  1 .  In  this  figure,  a  high  diffrac¬ 
tion  efficient  grating  is  used  at  plane  I*01  to  produce 
three  orders  of  smeared  color  spectra  at  the  back 
focal  plane  P„2-  Three  pinholes  (one  with  a  blue  filter) 
arc  properly  placed  over  this  set  of  spectral  lights,  to 
produce  three  primary  color  point  sources  (i.e.,  spat¬ 
ially  small)  at  plane  Pq2.  In  the  Fourier  transform 
plane  P2,  three  color  sensitive  inverse  filters  are  used 
for  the  restoration,  and  the  deblurred  color  images 
can  be  seen  at  the  output  plane  of  the  optical  proces¬ 
sor. 


2.  Dcfocuscd  color  image  restoration 

With  reference  to  the  optical  processor  of  fig,  I. 
there  arc  three  primary  color  point  sources  derived 
from  a  white-light  source.  The  intensity  distribution 
of  these  three  primary  color  point  sources  can  be 
represented  by  the  following  equation 

/(a.U;X)  =  A„cir|(cr2+/?2)'/2A/„| 

+  Xr  cir|((a  +  XR//>„)2  +  02)'/2/(/r  | 

+  Accir|((a-XG/p0)2+^)'/2H;|  .  (1) 

where  cir | X |  =  I ,  I  A|  <  I ,  and  0  otherwise;  (a.  / J )  is 
the  spatial  coordinate  system  of  plane  I’^j.p,,  is  the 
spatial  frequency  of  the  diffraction  grating,  T(x).f 
is  the  focal  length  of  the  achromatic  transform  lens, 
Ab,  Ar.Aq  are  the  proportional  constants,  </B,</R, 
d(~  arc  diameters  of  the  pinholes,  and  X„,  XR , 
are  the  center  wavelengths  of  the  blue,  red,  green  spec¬ 
tral  bands  of  the  light  sources,  respectively .  l-or  red 
and  green  point  sources,  the  center  wavelength  of 
each  spectral  baud  would  he 

X  =  a//P0.  (2) 

which  is  determined  by  the  position  of  the  pinhole 
at  a  axis.  The  corresponding  spectral  bandwidth  AX 
of  the  point  sources  can  be  written  as 
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Tig.  I .  A  white-lipht  optical  processor  for  out-of-focuscd  color  image  deblurring.  1;  wliite-liglu  source,  7'(jc);  diffraction  grating, 
R,  B.  G;  red,  blue,  and  preen  color  soutccs,  S(x,y)'.  out-of-focuscd  color  image,  L;  achromatic  transform  )enses,/(p,  Af;), 
ftp,  Ay), /(p,  Ar),  inverse  filters  lor  green,  blue,  red  light. 


AX  =  A Xt/fPn  =  (J/Jl>()  ,  (3) 

where  d  is  the  diameter  of  the  pinholes. 

We  note  that  the  complex  degree  of  coherence  pro¬ 
duced  by  the  color  point  sources  at  the  input  plane  Pj 
of  the  processor  is  113) 

2 f.(nd„r„) 

-  exp(i2n<£„x),  n  =  R,G,B  ,  (4) 

"unrn 

where  is  the  first-order  Bessel  function,  r„  = 
(l/A„/)(x2  +^2)'/2,and(/>„  =  0,-XRp0,and  XcP0 
for  the  blue,  red,  and  green  color  point  sources,  re¬ 
spectively.  Since  the  transform  lens  L02  is  achromatic, 
the  spatial  coherence  length  produced  by  a  circular 
point  source  would  be 

R=Kjyd,  (5) 

where  K  is  a  proportionality  constant,  and  K  -  1 .22 
if  the  spatial  coherence  length  R  is  defined  as  the  dis¬ 


tance  in  which  the  degree  of  coherence  p  drops  from 
unit  to  first  zero  value.  If  R  is  defined  from  a  drop  of 
12%  from  the  unit  degree  of  coherence,  then  K  -  0.32. 
With  reference  to  eqs.  (3)  and  (5),  we  note  that 

R  =  K\/AAp0  .  (6) 

Thus  to  maintain  the  same  spatial  coherence  length 
and  the  same  spectral  bandwidth,  different  spatial  fre¬ 
quency  Pq  of  the  grating  should  be  used.  However, 
with  the  use  of  multi-grating  frequency,  it  would  re¬ 
duce  to  available  power  of  the  light  source  for  the 
processing.  Nevertheless  this  problem  can  be  alleviated 
by  using  a  single  diffraction  grating  at  l*01  and  differ¬ 
ent  sizes  of  pinholes  at  ,,()2  plane.  In  this  manner  the 
same  spatial  coherence  length  can  be  obtained  by  sim¬ 
ply  varying  the  size  of  the  pinholes  of  the  primary 
color  sources.  As  a  numerical  example;  we  tabulate 
the  requirement  of  the  size  of  the  pinholes,  spatial 
coherence  lengths,  and  the  spectral  bandwidlhs  of  the 
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Table  I 

Relationship  or  si/.c,  spectral  bandwidth  and  spatial  coherence 
length  of  three  color  sources. 


X  (A) 

d  (pm) 

AX  (A) 

K 

R  (mm) 

6300 

400 

290 

0.32 

0.18 

red 

1.22 

0.67 

200 

140 

0.32 

0.35 

1.22 

1 .35 

5500 

400 

290 

0.32 

0.IS 

preen 

1.22 

0.59 

170 

120 

0.32 

0.36 

- 

1.22 

1.38 

4300 

400 

500 

0.32 

0.12 

blue 

1.22 

0.46 

140 

500 

0.32 

0.34 

1.22 

1.31 

Note:  AX  «  500  A  for  blue  wavelength  X  *  4300  A  is  limited 
by  Kodak  47B  filter. 


sources,  in  table  1 .  If  the  pinholes  for  red  and  preen 
sources  arc  of  the  saute  diameter,  e.g.,tyR  =  =  400 

pm,  the  spectral  bandwidths  would  be  AX  =  290  A. 

On  the  other  hand,  if  the  three  pinhole  diameters  arc 
different,  e.g.,t/R  =  200pm,</(;  =  1 70  pm,  and  t/},  = 
140  /mi, coherence  lengths  obtained  would  be  similar, 
i.e.,  1 .35  mm  (red),  1 .38  nun  (green),  and  1 .31  nun 
(blue). 

In  filler  synthesis,  we  mile  that  the  inverse  filler 
function  for  (lie  restoration  of  out-of-locus  image, 
for  each  primary  color  wavelength,  should  be 


/(P,A)  = 


27^)^ 

±  1  ,  otherwise  , 


where  a  is  the  diameter  of  the  defocused  point  spread 
function,  p  =  (l/X/)(a2  +  02)^2,  and  m  is  a  constant 
of  the  order  of  10  to  102.  The  meaning  of  m  is  that 
the  finite  filler  in  dynamic  range  can  be  made  in  phys¬ 
ical  sense  because  ./,(trop)  has  an  infinite  number  of 
zeros  Ml  .  Needless  to  say  (hat  cq.  (7)  can  also  be 
written  as 

f(p,  X)  =  l/(p,  X)|  c'°  . 

Tims  an  inverse  filter  can  be  synthesized  by  Ihc  com¬ 
bination  of  an  amplitude  and  a  phase  filter.  The  am¬ 
plitude  filter  can  be  synthesized  by  recording  an  am¬ 


plitude  spectrum  of  a  circular  aperture  of  a  given  diam¬ 
eter  using  a  lle-Nc  laser.  Since  the  scale  of  the  spec¬ 
trum  is  proportional  to  the  wavelength,  the  diameter 
of  the  circular  aperture  should  be  properly  fitted 
with  the  color  wavelength.  In  other  words,  the  diam¬ 
eter  for  Ihc  red  amplitude  filler  should  be  smaller 
than  the  green  amplitude  fillet,  and  the  diameter  lot 
the  green  filter  should  be  smaller  than  the  blue.  For 
example,  if  Ihc  diameter  of  the  circular  iqicrturc  for 
red  light  is  taken  equal  to  0.5  mm.  then  the  diameter 
of  Ihc  aperture  for  Ihc  green  should  be  0.58  mm.  and 
the  diameter  for  tlie  blue  tight  should  be  0.74  nun.  In 
amplitude  filter  synthesis,  it  is  necessary  to  control 
the  gamma  of  the  recording  plate  equal  to  1  to  obtain 
the  required  amplitude  transmittance. 

In  phase  filter  synthesis,  we  see  that,  the  filter  is 
primarily  composed  of  a  set  of  7r-phasc  circular  con¬ 
centric  rings,  as  shown  in  fig.  2.  The  diameters  of  the 
rings  are  determined  by  the  primary  color  wavelength 
and  the  si/.c  of  the  defocused  point  spread  function. 

A  simple  technique  of  producing  the  phase  filter  is 
a  bleaching  method  [14|.  A  black-and-white  ring  pat¬ 
tern,  corresponding  to  the  7r-phase  zone  of  the  phase 
filter,  is  recorded  on  a  high-contrast  film.  The  recorded 
binary  ring  pattern  is  used  as  a  mask  to  reproduce  a 
number  of  gray-level  rings  on  a  low-contrast  photo¬ 
graphic  plate.  If  the  recorded  plate  is  bleached,  a  set 
of  dcliliirring  phase  fillers  can  he  obtained.  To  search 
for  a  ff-phasc  filter  for  a  given  wavelength,  one  can 


-y,  n  phose  deloy 


*///A 


r„,  =  J?j.  Xfl! 
a 


where  )*o,  <wi,i,s,«, 

Xn-  XnISJnoU.  XjlSSOoi', 
X,i«  Jool), 

a  —  0.5  MM  .  t  —  550  MM  . 

Fig.  2.  A  n-phasc  concentric  rings  for  a  debturring  phase  filter. 
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utilize  (he  contrast-reverse  method  or  by  observing 
(he  intensity  ratio  of  Hie  main  peak  ami  lire  second 
peak  of  a  processed  point  spread  function. 

In  this  manner,  a  scl  of  demurring  niters  suitable 
for  the  three  primary  color  wavelengths  can  be  ob¬ 
tained.  The  primary  color  inverse  filter  functions  can 
be  described  as 


n  =  K.G.B , 


(8) 


where 

P„  4  (IAr/)U«  +  XR//»„)-  +  d’)'/2  .  for  XR  . 

(l/Xc/)Uft  -  V;/Po)2  +<32)'/2.  for  • 
(l/XB/)(<*2  +  02)!^  ,  for  XB  . 

The  three  inverse  filters  arc  placed  in  the  Courier  plane 
P2  of  the  processor,  their  positions  should  be  adjusted 
properly,  such  that  a  -  —  Xr/Po>  &  =  0  f°r  re^  sens*five 
Inverse  filler.  <«  ~  =  0  for  green,  and  tv  *  0, 

0  =  0  for  blue. 


3.  Experimental  results 

In  our  experiment,  a  75  W  xenon  arc  lamp  with  a 
500  /mi  pinhole  is  used  as  the  white-light  source  I. 

A  phase  giating  of  40  V/inm  Is  placed  at  plane  I',,) , 


and  three  pinholes  of  200  pm,  1 70 /am,  and  150 
pm,  for  the  red,  green,  and  blue  spectral  bands  are 
used  al  plane  •02  •  1  lie  center  wavelength  for  I  lie 
red  spectral  hand  is  6300  A,  5500  A  for  the  green, 
and  4300  A  for  the  blue.  Since  we  utilize  the  zero- 
order  spectra  for  the  blue  color  source,  a  Kodak  47B 
blue  color  filler  is  used  to  cover  the  pinhole.  The 
spccltal  bandwidtbs  of  these  ptiinaiy  color  turn  ices 
arc:  140  A  for  the  red,  1 20  A  for  the  green,  and 
500  A  Ini  the  trine.  Needless  to  say  Hint  Hie  Intensi¬ 
ties  of  these  three  primary  color  sources  can  be  ad¬ 
justed  by  a  set  of  neutral  density  filters  for  color 
balance. 

For  the  first  experimental  demonstration,  we  used 
a  circular  aperture  of  about  0.5  mm  in  diameter  as 
a  defocussed  point  spread  function,  as  shown  in  fig. 
3a.  Fig.  3b  shows  the  restored  image  of  fig.  3a  ob¬ 
tained  with  this  polychromatic  restoration  technique. 
From  this  figure,  we  sec  that,  the  dchlurrcd  point 
spread  function  mainly  consists  of  a  high  intensity 
center  peak,  as  the  dehlurred  image,  and  a  weak  cir¬ 
cular  ring  image,  which  is  primarily  due  to  the  finite 
extent  of  the  dcblurring  fillers.  In  our  experiment  the 
dcblurring  filters  arc  limited  to  about  four  lopes,  and 
the  greatest  density  of  the  amplitude  filter  is  about 
2  to  2.5  D. 

From  the  result  of  fig.  3b,  we  also  see  that  there 
is  a  slight  color  dispersion  at  two  edges  of  the  center 
restored  peak  image.  We  note  that,  the  dispeision  is 


t  ip.  3.  Restoration  of  a  dcfocused  point  spread  function,  a)  A  transparent  circular  disk  as  a  dcfocuscd  point  spread  function, 
b)  A  black-and-white  picture  of  the  dcblurrcd  point  spread-function. 
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Fig.  4,  A  electronic  scanned  image  of  fig.  3b. 

primarily  due  to  the  chromatic  aberration  of  the 
transform  lenses,  which  can  be  eliminated  by  using 
higher  (|uulity  achromatic  transform  lenses.  Fig.  4 
shows  an  electronic  scanned  image  of  the  result  of 
fig.  3b.  From  this  Figure,  we  see  that  the  intensity  of 
the  center  peak  is  much  higher  than  that  of  the  first- 
order  ring. 

We  now  provide  a  second  experimental  demonstra¬ 
tion  for  three  overlapping  primary  color  disks,  as 
shown  in  fig.  5a.  The  size  of  these  primary  color  disks 


Fig.  5.  Restoration  of  three  defocused  color  point  spread  func¬ 
tions.  a)  A  black-and-white  photograph  of  three  overlapping 
primary  color  disks  as  the  color  defocused  point  spread  func¬ 
tions.  b)  A  black-and-white  photograph  of  the  corresponding 
dcblurred  point  spread  functions. 


arc  about  0.5  mm  in  diameter.  If  the  color  transparen¬ 
cy  of  Fig.  5a  is  inserted  at  the  input  plane  of  the  pro¬ 
posed  deblurring  processor  of  Fig.  I ,  the  restored  color 
point  images  can  be  obtained,  as  shown  in  Fig.  5b. 

From  this  Figure  we  sec  that  the  color  of  the  restored 
images  are  very  faithful.  However  there  is  a  slight  re¬ 
placement  of  the  dcblurred  color  point  image,  which 
is  primarily  due  to  the  chromatic  aberration  of  the 
transform  lenses. 

As  a  final  experimental  result,  fig.  6a  shows  two 
blurred  color  words  (i.e.,  “Color  Image”)  as  an  input 
binned  image.  The  cnnespnitding  debit! tied  color 
image  obtained  with  this  dchhirring  lcclmi(|uc  is  shown 
in  fig.  6b.  From  this  deblurred  color  image,  again  we 
see  that  the  color  reproduction  is  rather  faithful  and 
the  quality  of  the  dcblurred  image  is  quite  impressive. 
Since  the  color  image  dcblurring  is  obtained  with  three 


c 
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Fig.  6.  Color  image  restoration  of  defocused  color  words 
"Color  Image”,  a)  A  black-and-white  photograph  of  the  de¬ 
focused  color  words,  b)  A  black-and-white  photograph  of  the 
corresponding  dcblurred  color  image. 
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relatively  temporally  broad  and  spatially  large  prima¬ 
ry  color  light  sources,  the  finger  print  like  structure 
coherent  artifact  noise  is  eliminated,  as  can  be  seen 
from  figs.  3b,  5b,  and  6b. 


4.  Summary 

We  have  shown  that  the  color  image  restoration  for 
out-of-focuscd  photographic  images  can  be  obtained 
by  a  white-light  processing  technique.  This  technique 
uses  a  diffraction  grating  to  generate  three  diffraction 
orders,  so  that  three  relatively  spectrally  broad  and 
spatially  large  primary  color  sources  can  be  derived 
from  these  diffractions.  Since  the  temporal  and  spat¬ 
ial  coherence  requirements  for  the  image  dcblurring 
are  determined  by  the  spectral  bandwidth  and  the 
spatial  size  of  the  primary  color  sources,  the  three 
primary  color  partially  coherent  sources  can  be  ob¬ 
tained  by  this  proposed  technique.  The  advantage  of 
the  technique  is  that  one  can  adjust  the  degree  of  the 
temporal  and  spatial  coherence  of  each  primary  color 
source  by  simply  changing  its  diameter.  By  utilizing 
these  primary  color  sources,  an  out-of-focused  blurred 
color  photographic  image  can  be  restored  by  a  set  of 
color  sensitive  inveisc  fillers.  And  good  dcblurred  color 
images  have  been  obtained  with  pinhole  sizes  as  large 
as  400  /am.  To  alleviate  the  low  diffraction  efficiency 
of  the  holographic  inverse  fillers,  the  filter  synthesis 
arc  obtained  by  the  combination  of  non-absorptive 
phase  filters  and  absorptive  amplitude  filters.  The  phase 
filters  are  obtained  with  a  bleaching  technique,  while 
the  amplitude  filters  are  obtained  by  intensity  expo¬ 
sure. 


In  view  of  the  experimental  results,  we  see  that  the 
resolved  color  images  offer  a  reasonably  good  deblurred 
image  quality  and  the  colors  are  faithfully  reproduced. 
Since  the  primary  color  light  sources  are  spectrally  and 
spatially  broad,  the  coherent  artifact  noise  is  substan¬ 
tially  suppressed.  Although  there  is  some  degree  of 
chromatic  aberration,  it  can  be  alleviated  by  utilizing 
higher  quality  achromatic  transform  lenses. 
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SOURCE  ENCODING,  SIGNAL  SAMPLING 
AND  SPECTRAL  BAND  FILTERING 
FOR  PARTIALLY  COHERENT  OPTICAL  SIGNAL 

PROCESSING 
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Mots  cl£s 

Traitement  optique 
Coherence 


Traitement  optique  en  cclairage  particllcmcnt  coherent 
en  modulant  la  source  et  le  signal 
et  en  filtrant  le  spectre  du  signal 

RfeSUMfe  :  On  disculc  dcs  relations  cnirc  la  coherence,  le  mirage. 
I'cchantillonnagc  du  signal  el  la  modulation  dc  la  source.  Puisquc 
le  degre  dc  coherence  spalialc  cst  determine  par  le  traitement  dc 
Fintage.  tine  parfailc  coherence  spalialc  n'est.  en  general,  pas  ncecs- 
SiiirC.  L'avantagc  dc  modulcr  la  source  cst  de  permctlrc  dc  trailer 
f  image  avee  unc  source  6lenduc  incohercnlc.  L'cffcl  d’cchunlil- 
lonner  I'image  est  dam6liorcr  la  coherence  tcmporelle  dans  le  plan 
de  Fourier  de  facon  que  le  filtrage  puisse  s'elTectuer  en  lumiere 
parliellement  coherente  temporellement.  L' utilisation  d'un  domaine 
spectral  large  a  pour  but  dc  rtduire  le  bruit  du  a  la  coherence 
spatiale  de  I'cclairagc.  F.nlin.  puisque  un  disposilif  dc  traitement 
d'imagc  en  lumiere  particllcmcnt  coherente  utilise  unc  source  dc 
lumiere  blanche,  il  purait  particuliercmcnl  bicn  adapte  au  traite¬ 
ment  des  images  en  couleurs.  On  presente  des  experiences  mon- 
trant  l'avantagc  de  modulcr  la  source  ct  d’£chantillonncr  I'image 
lorsqu'on  utilise  unc  source  ctcnduc  dc  lumiere  blanche. 


SUMMARY  :  Relations  between  coherence  requirement,  spectral 
filtering  signal  sampling  and  source  encoding  arc  discussed.  Since 
the  spatial  coherence  requirement  is  determined  by  the  signal  pro¬ 
cessing  operation,  a  strict  spatial  coherence  is  usually  not  required. 
The  advantage  of  the  source  encoding  is  to  relax  the  constraints 
of  a  physical  light  source  so  that  the  signal  processing  can  be  carried 
out  with  an  extended  incoherent  source.  The  effect  of  signal  sampling 
is  to  improve  the  temporal  coherence  requirement  at  the  Fourier 
plane  so  that  the  spatial  filtering  can  be  carried  out  with  partially 
coherence  mode.  The  objective  of  broad  spectral  band  filtering  is 
to  carry  out  the  signal  processing  over  the  entire  spectral  band  of 
the  light  source  so  that  the  coherent  noise  can  be  eliminated.  Since 
the  partially  coherent  optical  processor  utilizes  a  broad  spectral 
band  white-light  source,  it  is  particularly  suitable  for  color  signal 
processing.  Experimental  demonstrations  for  the  source  encoding 
signal  sampling  and  spectral  band  filtering  arc  included. 


Key  words 

Optical  processing 
Coherence 


INTRODUCTION 

Since  the  invention  of  laser  (i.c..  u  strong  coherent 
source)  laser  has  become  a  fashionable  tool  for  many 
scientific  applications  particularly  as  applied  to  cohe¬ 
rent  optical  signal  processing.  However  coherent 
optical  signal  processing  systems  are  plagued  with 
coherent  noises,  which  frequently  limit  their  process¬ 
ing  capability.  As  noted  by  the  late  Gabor,  the  Nobel 
prize  winner  in  physics  in  1970  for  his  invention  of 
holography,  the  coherent  noise  is  the  number  one 
enemy  of  the  Modern  Optical  Signal  Processing  (1). 
Aside  the  coherent  noise,  the  coherent  sources  arc 
usually  expensive,  and  the  coherent  processing  envi¬ 


ronments  arc  very  stringent.  For  example,  heavy 
optical  benches  and  dust  free  environments  are  gene¬ 
rally  required. 

Recently,  we  have  looked  at  the  or.«:al  processing 
from  a  different  standpoint.  A  qu.c.v.i  arises,  is  it 
necessarily  true  that  all  optical  signal  processing 
required  a  coherent  source  ?  The  answer  to  this 
question  is  that  there  arc  many  optical  signal  pro¬ 
cessings  that  can  be  carried  out  by  a  white-light 
source  [2].  The  advantages  of  the  proposed  white- 
light  signal  processing  technique  are  :  1.  It  is  capable 
of  suppressing  the  coherent  noise;  2.  White-light 
sources  arc  usually  inexpensive;  3.  The  processing 
environments  arc  not  critical;  4.  The- white-light  sys- 
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tcm  is  relatively  easy  and  economical  to  maintain; 
and  5.  The  white-light  processor  is  particularly  sui¬ 
table  for  color  image  processing.  - 

One  question  that  the  reader  may  ask.  since  the 
white-light  system  offers  all  these  glamorous  merits, 
why  it  has  been  ignored  for  so  long  ?  The  answer  to 
this  question  is  that,  it  was  a  general  acceptance  that 
an  incoherent  source  cannot  process  the  signal  in 
complex  amplitude.  However,  none  of  the  practical 
sources  are  strictly  incoherent,  even  a  white-light 
source.  In  fact,  we  were  able  to  utilize  the  partial 
coherence  of  a  white-light  source  to  perform  the 
complex  amplitude  processing.  The  proposed  white- 
light  processor,  on  one  hand  it  is  capable  of  suppressing 
the  coherent  noise  like  an  incoherent  processor,  on 
the  other  hand  it  is  capable  of  processing  the  signal 
in  complex  amplitude  like  a  coherent  processor. 

There  is  however  a  basic  different  approach  toward 
a  coherent  and  a  white-light  processor.  In  coherent 
processing,  virtually  no  one  seems  to  care  about  the 
coherence  requirements,  since  the  laser  provides  a 
strong  coherent  source.  However,  in  white-light  pro¬ 
cessing.  the  knowledge  of  the  coherence  requirement 
is  usually  needed. 

In  white-light  processing  we  would  approach  the 
problem  backward.  First,  we  should  know  what  is  the 
processing  operation  we  wish  to  perform  :  Is  it  a  I -D 
or  2-D  processing  ?  Is  the  signal  filtering  a  point  or 
point-pair  concept  ?  What  is  the  spatial  bandwidth 
of  the  signal  ?  etc.  Then  with  these  knowledges,  we 
would  be  able  to  evaluate  the  coherence  requirements 
at  the  Fourier  and  at  the  input  planes.  From  the  eva¬ 
luated  results,  we  would  be  able  to  design  a  signal 
sampling  function  and  a  source  encoding  function 
to  obtain  these  requirements.  The  objective  of  using 
a  signal  sampling  function  is  to  achieve  a  high  degree 
of  temporal  coherence  in  Fourier  plane  so  that  the 
signal  can  be  processing  in  complex  amplitude,  for 
the  entire  spectral  band  of  a  white-light  source.  And 
for  the  source  encoding  is  to  alleviate  the  constraint 
of  an  extended  white-light  source. 

In  the  following  sections,  we  shaH  discuss  in  detail 
the  source  encoding,  signal  sampling  and  spatial  band 
filtering  as  applied  to  a  partially  coherent  optical 
(c.g.,  white-light)  signal  processing. 

PARTIALLY  COHERENT  OPTICAL  SIGNAL 
PROCESSING 

Wc  shall  now  describe  an  optical  signal  processing 
technique  that  can  be  carried  out  by  a  broad  band 
white-light  source,  as  illustrated  in  figure  /.  The 
white-light  signal  processing  system  is  similar  to  that 
of  a  coherent  system,  except  the  use  of  a  white-light 
source,  source  encoding  mask,  signal  sampling  grating, 
multispcctral  filters  and  achromatic  transform  lenses. 
For  example,  if  wc  place  a  signal  transparency  s(x,  y) 
in  contact  with  a  sampling  phase  grating,  the  complex 
light  field  for  every  wavelength  A  behind  the  achro¬ 
matic  transform  lens  L,  would  be 

E(p.  q ;  A)  =  J  J  rfx.  y)  exp(i/>0  x)  x 

x  cxp[-  i(/).v  +  </»•)]  d.v  dr  =  Sip  -  p0,  </) ,  (1) 


where  the  integral  is  over  the  spatial  domain  of  the 
input  plane  /’,.  (/>,</)  denotes  the  angular  spatial 
frequency  coordinate  system,  p0  is  the  angular  spatial 
frequency  of  the  sampling  phase  grating,  and  S(p.  q) 
is  the  Fourier  spectrum  of  x(x.y).  If  wc  write  F.q.  (I) 
in  the  form  of  linear  spatial  coordinate  system  (a,  /I), 
we  have 

E(u,p-,X)  =  s(^-jlPoJty  (2) 

where  p  -  (2  n/X.f)  a,  q  ^  (2  n/X.f)  fi  and  /  is  the 
focal  length  of  the  achromatic  transform  lens.  Thus, 
we  see  that  the  Fourier  spectra  would  disperse  into 
rainbow  color  along  the  a  axis,  and  each  Fourier 
spectrum  for  a  given  wavelength  X  is  centered  at 
a  ±  (XJJ2  n)  p0. 


Extended  » 
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Fio.  I .  —  A  whiie-ligltl  optical  signal  processor. 


In  signal  filtering,  wc  assume  that  a  sequence  of 
complex  spatial  filters  for  various  A*  are  available. 
i.e.,  //(/>„.  qn).  where  pn  =(2  n/XJ)  a,  q„  =(2  njXJ)  0. 
In  practice,  all  the  processing  signals  arc  spatial  fre¬ 
quency  limited,  the  spatial  bandwidth  or each  spectral 
band  filter  //(/»„,  q„)  is  also  bandlimitcd.  such  as 


„>  <?»)  =  | 


H(p„qn),  a,  <  a  <  a2 , 
0 ,  otherwise , 


where  a,  ^  (A„/72  n)(p0  +  A p)  and  oc2  =s(XKfi2  n)  x 
( p0  -  A p)  are  the  upper  and  the  lower  spatial  limits 
of  //(p„,  </„).  and  A p  is  the  spatial  bandwidth  of  the 
input  signal  ,v(x,  y ). 

The  limiting  wavelengths  of  each  //(/>„,  </„)  can  be 
written  as 

,  ,  /’»  +  A/I  „  ,  ,  ,  Pn  -  A/>  IAI 

X.  =  A_ - j-.  and  -  X„ - — r— .  (4) 

'  "  P„  -  A p  />„  i  A/ 1 

The  spectral  bandwidth  of  //(/v  q„)  is  therefore. 

H-A.  V  (5) 

’  P1  ~  (A p)1  Pn 

If  we  place  this  set  of  spectral  band  filters  side-by- 
side  positioned  over  the  smeared  Fourier  spectra,  then 
the  intensity  distribution  of  the  output  light  field  can 
bi  shown  as. 

fix.  y)  *  £  A  A.  |  s(x,  y;  A„)  *  fi(x,  y;  A„)  |2 ,  (6) 

nr-  I 

where  /i(.v,  y ;  A)  is  the  spatial  impulse  response  of 
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ll[pn.  </„)  and  *  denotes  the  convolution  operation. 
Thus,  the  proposed  while-light  signal  processor  is 
capable  of  processing  the  signal  in  complex  amplitude. 
Since  the  output  intensity  is  the  sum  of  the  mutually 
incoherent  narrow  hand  irradianees.  the  annoying 
coherent  noise  can  be  eliminated.  Furthermore,  the 
white-light  source  contains  all  the  color  wavelengths, 
the  proposed  system  is  particularly  suitable  for  color 
signal  processing. 

SPECTRAL  BAND  FILTERING,  SIGNAL 
SAMPLING  AND  SOURCE  ENCODING 

We  have  mentioned  earlier  for  white-light  or  par¬ 
tially  coherent  processing,  we  would  approach  the 
problem  in  backward  manner.  For  example,  if  signal 
filtering  is  two-dimensional  (c.g;.  2-1)  correlation  ope¬ 
ration),  we  would  synthesize  a  set  of  narrow  spectral 
band  filters  for  each  for  the  entire  smeared  Fourier 
spectra,  as  illustrated  in  figure  2(a).  On  the  other 
hand,  if  the  signal  filtering  is  one-dimensional  (e.g, 
deblurring  due  to  linear  motion),  a  broadband  fan- 
shape  spatial  filter,  to  accomodate  the  scale  variation 
due  to  wavelength,  can  be  utilized  as  illustrated  in 
figure  2(b).  Since  the  filtering  is  taken  place  with  the 


Fk;  2  —  (a)  A  multi  .ipertrul-huml Jiltrr.  (b)  A  fun-sluipr  filler. 
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entire  spectral  band  of  the  light  source,  the  coherent 
noise  can  be  suppressed  and  the  white-light  processing 
technique  is  very  suitable  for  colour  image  processing 
There  is.  however,  a  temporal  coherence  require¬ 
ment  imposed  upon  the  signal  filtering  in  Fourier 
plane.  Since  the  scale  of  the  Fourier  spectrum  varies 
with  wavelength,  a  temporal  coherence  requirement 
should  be  imposed  on  each  spatial  filter  at  the  Fourier 
plane.  Thus,  the  spectral  spread  over  each  filter 
H(p„.  qn)  is  imposed  by  the  temporal  coherence  requi¬ 
rement,  i.e.. 


(7) 


From  this  requirement,  a  high  degree  of  temporal 
coherence  is  achievable  in  the  Fourier  plane  by  simply 
increasing  the  spatial  frequency  of  the  sampling  grat¬ 
ing.  Needless  to  say  that  the  same  temporal  coherence 
requirement  of  Eq.  (7)  can  also  be  applied  for  a 
broadband  fan-shape  filter. 

There  is  also  a  spatial  coherence  requirement 
imposed  at  the  input  plane  of  the  white-light  processor. 
With  reference  to  the  Wolfs  [3]  partial  coherence 
theory  |3).  the  spatial  coherence  function  at  the  input 
plane  can  be  shown  [4], 


f(x  -  x')  =  j*j"  y(x0)  expj^i  2  n  ^(x  -  x')  dx0 


where  y(x0)  denotes  the  intensity  distribution  of  the 
source  encoding  function. 

From  the  above  equation,  wc  see  that  the  spatial 
coherence  and  source  encoding  functions  form  a 
Fourier  transform  pair,  i.e.. 


y(x0)  =  .T[r(x  -  x’)].  (9) 

and 


r(\-x’)  =  j  1  [y(x0)] ,  (io) 


where  .7  denotes  the  Fourier  transformation.  This 
Fourier  transform  pair  implies  that  if  a  spatial  cohe¬ 
rence  function  is  given  then  the  source  encoding 
function  can  be  evaluated  through  the  Fourier  trans¬ 
formation  and  vice  versa.  Wc  note  that  source  encoding 
function  can  consist  of  apertures  of  any  shape  or 
complicated  gray  scale  transparency.  However  the 
source  encoding  function  is  only  limited  to  a  positive 
real  quantity  which  is  restricted  by  the  following 
physical  realizable  condition  : 


0  <  y(x„)  I  .  (II) 

In  white-light  processing,  wc  would  search  for  a 
reduced  spatial  coherence  requirement  for  the  pro¬ 
cessing  operation.  With  reference  to  this  reduced 
spatial  coherence  function,  a  source  encoding  func¬ 
tion  that  satisfied  the  physical  realizability  condition 
can  be  obtained.  One  of  the  basic  objectives  of  the 
source  encoding  is  to  alleviate  the  constraint  of  a 
white-light  source.  Furthermore  the  source  encoding 
also  improves  the  utilization  of  the  light  power  such 
that  the  optical  processing  can  be  carried  out  by  an 
extended  source. 
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Ho.  3.  —  A  wlute-liglii  image  subtraction  processor,  7(.v);  phase 
grating.  L, ;  image  lens.  l.f  ,  collimated  lens,  Lt  and  L, ;  achromatic 
iriiiis/ortH  lenses,  y(  i) .  source  entailing  mask,  (i .  hm-shape  dij- 
J taction  grating. 


We  shall  now  illustrate  an  application  of  the  source 
encoding,  signal  sampling  and  Tillering  for  a  white- 
light  signal  processing.  Let  us  now  consider  a  poly¬ 
chromatic  image  subtraction  |5).  The  image  sub¬ 
traction  of  Lee  [ft]  that  we  would  consider  is  essentially 
a  one-dimensional  processing  operation,  in  which  a 
1  -D  fan-shape  diffraction  grating  should  be  utilized, 
as  illustrated  in  figure  3.  We  note  that  the  fan-shape 
grating  (i.e..  filter)  is  imposed  by  the  temporal  cohe¬ 
rence  condition  of  Eq.  (7).  Since  the  image  subtraction 
is  a  point-pair  processing  operation,  a  strictly  broad 
spatial  coherence  function  at  the  input  plane  is  not 
required.  In  other  words,  if  one  maintains  the  spatial 
coherence  between  the  corresponding  image  points 
to  be  subtracted  at  the  input  plane,  then  the  subtrac¬ 
tion  operation  can  be  carried  out  at  the  output  image 
plane.  Thus  instead  of  using  a  strictly  broad  spatial 
coherence  function,  a  reduced  spatial  coherence  func¬ 
tion  may  be  utilized,  such  as 

r(.v  -  v')  =  tff.v  -  .»•'  -  /in)  +  S(y  -  y  +  h0),  (12) 

where  2 /i0  is  the  main  separation  between  the  two 
input  color  transparencies.  The  source  encoding  func¬ 
tion  can  therefore  be  evaluated  by  through  the  Fourier 
transform  of  Eq.  (9).  such  as 

y(  i  0)  =  2  cos  •  (13) 

Unfortunately  Eq.  (13)  is  a  bipolar  function  which  is 
not  physically  realizable.  To  ensure  a  physically 
realizable  source  encoding  function,  we  let  a  reduced 
spatial  coherence  function  with  the  required  point- 
pair  coherence  characteristic  be  (7]. 


/  (I  >’-/!)  = 
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where  NH  a  positive  integer,  and  it  4  <1.  Eq.  (14) 
represents  a  sequence  of  narrow  pulses  which  occur 
at  every  |  y  -  y' 1  =  «/r0,  where  n  is  a  positive  integer, 
and  their  peak  values  arc  weighted  by  a  broader  sine 
factor,  as  shown  in  figure  4(a).  Thus,  a  high  degree 
of  spatial  coherence  can  be  achieved  at  every  point- 
pair  between  the  two  input  color  transparencies.  By 
taking  the  Fourier  transformation  of  the  reduced 
spatial  coherence  function  of  Eq.  (14),  the  correspond¬ 
ing  source  encoding  function  is 

7(1  J’D-  I  rcct  nd\  (|5) 

■  -  I  W 

where  w  is  the  slit  width,  d  =  (A///»0)  is  the  separation 
between  the  slits,  and  N  is  the  number  of  the  slits. 
Since  y(|  y  f)  is  a  positive  real  function  which  satisfies 
the  constraint  of  Eq.  (II),  the  proposed  source  encod¬ 
ing  function  of  Eq.  (15)  is  physically  realizable. 

In  view  of  Eq.  (1 5)  we  also  note  that,  the  separation 
of  slit  d  is  linearly  proportional  of  the  A.  The  source 
encoding  is  a  fan-shape  type  function,  as  shown  in 
figure  4(h).  To  obtain  lines  of  rainbow  color  spectral 
light  source  for  the  signal  processing,  we  would  utilize 
a  linear  extended  white-light  source  with  a  dispersive 
phase  grating,  as  illustrated  in  figure  3.  Thus  with 


r<»-») 
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Flo.  4.  —  (n)  A  spatial  coherent  c  turn  lion  (h)  4  source  mauling 
mask . 
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appropriate  somvc  encoding,  signal  dimpling  and  fil- 
Iciing.  color  image  suhliaclion  opcialion  can  Iv 
obtained  at  the  output  plane.  We  stress  again,  the 
basic  advantage  of  source  encoding  is  to  alleviate  the 
constraint  of  strict  spatial  coherence  requirement 
imposed  upon  the  optical  signal  processor.  The  source 
encoding  also  offers  the  advantage  of  efficient  utili¬ 
zation  of  the  light  power. 


EXPERIMENTAL  DEMONSTRATIONS 

We  shall  now  provide  a  couple  of  experimental 
results  obtained  with  the  source  encoding,  signal 
sampling  and  spectral  band  filtering  technique  for 
while-light  and  extended  incoherent  courses.  We 
shall  first  show  the  result  obtained  for  color  image 
dcblurring  due  to  linear  motion  with  the  white- 
light  processing  technique.  Since  linear  motion  deblur¬ 
ring  is  a  l-D  processing  operation  and  the  inverse 
filtering  is  a  poinl-by-poinl  filtering  concept  such 
that  the  operation  is  taking  place  on  the  smearing 
length  of  the  blurred  image.  Thus  the  deblurring 
filter  (i.c..  inverse  filler)  is  a  fan-shape  type  spatial 
filler  | SJ  ami  the  temporal  coherence  requirement 
is  imposed  by  Fq.  (7).  1  he  spatial  coherence  require¬ 
ment  is  dependent  upon  the  smearing  length.  A  source 
encoding  function  of  a  narrow  slit  width  (dependent 
upon  the  smearing)  perpendicular  to  the  smearing 
length  is  utilized.  Figure  5(a)  shows  a  color  picture 
of  a  blurred  image  due  to  linear  motion  of  a  F-16 
fighter  plane.  The  body  of  this  fighter  plane  is  painted 
in  navy  hluc-und-whitc  colors,  the  wings  are  mostly 
painted  in  red.  the  tail  is  navy  bluc-and-whilc.  and  the 
ground  terrain  is  generally  bluish  color.  From  this 
figure,  we  see  that  the  plane  is  badly  blurred. 
Figure  5(h)  shows  the  color  image  dcblurring  result 
that  we  obtained  with  the  proposed  while-light 
dcblurring  technique.  From  this  deblurrcd  result, 
the  letters  and  overall  shape  of  the  entire  airplane 
are  more  distinctive  than  the  blurred  one.  Further¬ 
more  the  river,  the  highways,  and  the  forestry  of  the 
ground  terrain  arc  far  more  visible.  We  note  that 
the  color  reproduction  of  the  deblurrcd  image  is 
spectacularly  faithful,  and  coherent  artifact  noise  is 
virtually  non-cxistcd.  There  is.  however,  some  degree 
of  color  blur  and  color  deviation,  w  hich  arc  primarily 
due  to  the  chromatic  aberration  and  the  anti-reflec- 
tancc  coaling  of  the  transform  lenses.  Neverthe¬ 
less.  these  drawbacks  can  be  overcome  by  utilizing 
good  quality  nchtomntic  tiansform  lenses. 

let  us  now  provide  a  color  image  subtraction 
utilized  by  the  source  encoding  technique  with 
extended  incoherent  sources  as  described  in  previous 
sections.  Figure  6(a)  and  6(h)  show  two  color 
image  transparencies  of  a  parking  lot  as  input  color 
objects.  Figure  6(e)  shows  the  color  subtracted 
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image  obtained  by  the  source  encoding  technique 
with  extended  incoherent  source.  In  this  figure, 
the  profile  of  a  (red)  subcompacl  car  can  be  seen  at 
the  output  image  plane.  The  shadow  and  the  parking 
line  (in  yellow  color)  can  also  be  readily  identified. 
We  however  note  that,  this  color  image  subtraction 
result  is  obtained  by  two  narrow  band  extended 
incoherent  sources.  Extention  toward  the  entire 
spectral  band  of  a  white-light  source  is  currently 
under  investigation. 

CONCLUSION 

In  conclusion  we  would  point  out  that  the  advantage 
of  source  encoding  is  to  provide  an  appropriate 
spatial  coherence  function  at  the  input  plane  so  that 
the  signal  processing  can  be  carried  out  by  an  extended 
incoherent  source.  The  effect  of  the  signal  sampling 
is  to  achieve  the  temporal  coherence  requirement 
at  the  Fourier  plane  so  that  the  signal  can  be  processed 
in  complex  amplitude.  If  the  filtering  operation  is 
two-dimensional,  a  multi-spectral-band  2-1)  filters 
should  be  utilized.  If  the  filtering  operation  is  one- 
dimensional.  a  fan-shape  filter  can  be  used. 

In  short,  one  should  carry  out  the  processing 
requirements  backward  for  a  partially  coherent  or 
white-light  processing.  With  these  processing  require¬ 
ments  (eg.,  operation,  temporal  and  spatial  coherence 
requirements),  multi-spcctral-band  or  fan-shape  filter, 
signal  sampling  function,  and  source  encoding  mask 
can  be  synthesized.  Thus  the  signal  processing  can 
be  carried  out  in  complex  amplitude  over  the  whole- 
spectral  band  or  an  extended  white-light  source. 

We  acknowledge  the  support  of  the  U.S.  Air  Force 
Office  of  Scientific  Research  Grant  AFOSR-8I-OI48. 
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Advances  in  White-Light  Signal  Processing 
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SUMMARY 

A  technique  Lhnt  permits  signal  processing  operations  to  be  carried  out  by 
white-light  source  is  described.  This  method  is  capable  of  performing  signal, proc¬ 
essing  that  obeys  the  concept  of  coherent  light  ratiirr  than  incoherent  option.  Since 
the  white- .light  source  contains  all  the  color  wavelengths  of  the  visible  light,  the 
technique  is  very  suitable  for  color  signal  processing.  Some  of  its  recent  advances 
in  white-light  signal  processing  will  be  illustrated. 

INTRODUCTION 

The  use  of  coherent  light  enables  optical  systems  to  carry  out  many  sophisti¬ 
cated  information  processing,  operations  (ref.  1).  However,  coherent  optical  proc¬ 
essing  systems  are  plagued  with  coherent  artifact  noise,  which  freoucntly  limits 
their  processing  capability.  Although  many  optical  information  processing  operations 
can  be  implemented  by  systems  that  use  incoherent  light  (refs.  2-5),  there  are  other 
severe  drawbacks.  The  incoherent  processing  system  Is  capable  of  reducing  the  inev¬ 
itable  artifact  noise,  but  it  generally  introduces  a  dc-bias  buildup  problem,  which 
results  in  poor  noise  performance.  Techniques  have  been  developed  for  coherent 
operation  with  light  of  reduced  coherence  (refs.  6,7);  however,  these  teebnioues 
also  possess  severe  limitations. 

Attempts  at  reducing  the  temporal  coherence  reouirements  on  the  light  source  in 
optical  information  processing  fall  into  two  general  categories:  one,  the  use  of 
incoherent  instead  of  coherent  optical  processing  has  been  pursued  by  Loventhal  and 
Chavcl  (ref.  8)  and  Lolimann  (ref.  9),  among  others.  The  other,  the  reduction  of 
coherence  while  still  operating  in  the  linear-in-amplitude,  has  been  pursued  by 
Leith  and  Roth  (ref.  10)  and  by  Morris  and  George  (ref.  11). 

Since  the  invention  of  laser  (i.e.,  a  strong  coherent  source),  it  has  become 
a  fashionable  tool  for  many  scientific  applications  particularly  ns  applied  to 
coherent  optical  signal  processing.  However  coherent  optical  signal  processing 
systems  arc  plagued  with  coherent  noises,  which  frequently  limit  their  processing 
capability.  As  noted  by  the  late  Gabor,  the  Nobel  prize  winner  in  physics  in  1970 
for  his  invention  of  holography,  the  coherent  noise  is  the  number  one  enemy  of  the 
Modern  Optical.  Signal  1’roceMii  lug  (ref.  12).  Aside  Crum  tile  coherent,  no  too,  the 
coherent  Hourcrn  are  usually  expensive,  and  the  coherent  processing  environments 
'  «re  very  stringent.  For  example,  heavy  optical  benches  and  dust  free  environments 
are  generally  required. 

Recently,  we  have  looked  at  the  optical  processing  from  a  different  standpoint. 

A  question  arises,  is  it  necessarily  true  that  all  optical  signal  processing  required 

•This  research  is  supported  in  part  by  the  Air  Force  Office  of  Scientific  Research. 
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;i  cnluTiMtt  fiimri’f?  Tin;  answer  Ln  tliln  question  Jm  tli.it  l  here  are  many  optical 
signal  process  lugs  tlint  ran  be  carried  out  by  a  while-light  source  (ref.  IT).  The 
advantages  of  the  proposed  whlle-J  iglil  signal.  processing  lechnlaue  are:  1.  It  Is 
capable  of  suppressing  the  coherent  noise;  2.  White-light  source  is  usually 
inexpensive;  3.  The  processing  environment  is  not  very  demanding;  4.  The  white- 
light  system  is  relatively  easy  and  economical  to  maintain;  and  5.  The  white-light 
processor  is  particularly  suitable  for  color  image  processing. 

One  question  that  the  reader  may  ask,  since  the  white-light  system  offers  all 
these  glamorous  merits,  is  why  has  it  been  ignored  for  so  long?  The  answer  to  this 
question  is  that  it  was  a  general  acceptance  that  an  incoherent  source  cannot 
process  the  signal  in  complex  amplitude.  However,  none  of  the  practical  sources  are 
strictly  incoherent,  even  a  white-light  source.  Tn  fact,  we  were  able  to  utilize 
the  partial  coherence  of  a  wh  I  t«*-l  Jglit  source  Lo  perform  the  complex  amplitude 
processing.  The  proposed  white-light  processor,  on  one  hand,  is  capable  of 
suppressing  the  coherent  noise  like  an  Incoherent' .processor,  and  on  the  other  hand.  It 
is  capable  of  .processing  the  signal  in  complex  amplitude  like  a  coherent  processor. 

There  is  however  a  different  approach  toward  the  utilization  of  a  white-light 
processor.  In  coherent  processing,  virtually  no  one  evaluates  the  coherence 
requirement,  since  the  laser  provides  a  very  good  coherent  source.  However,  in 
white-light  processing,  the  evaluation  of  the  coherence  requirement  is  usually 
called  for. 

In  white-light  signal  processing  we  should  annroaeb  the  oroblem  from  a  different 
standpoint.  First,  we  should  have  the  a  priori  knowledge  of  the  signal  processing 
operation  we  would  encounter.  For  example,  is  it  a  1-D  or  2-D  processing?  Is  the 
signal  filtering  a  point  or  point-pair  concept?  What  is  the  spatial  bandwidth  of 
the  signal?  Then  we  would  be  able  to  evaluate  the  coherence  requirements  at 
the  Fourier  and  at  the  input  planes.  From  the  evaluated  results,  we  would  be  able 
to  design  a  signal  sampling  function  and  a  source  encoding  mask  to  obtain  these 
requirements.  The  objective  of  using  a  signal  sampling  function  is  to  achieve  a 
high  degree  of  temporal  coherence  in  Fourier  plane  so  that  the  signal  can  be  proc¬ 
essed  in  complex  amplitude  for  the  entire  spectral  band  of  the  light  source.  And 
for  the  source  encoding,  it  is  to  alleviate  the  inability  of  an  extended  source. 

In  the  following,  we  shall  discuss  in  detail  the  source  encoding,  signal 
sampling  and  spatial  band  filtering  as  applied  to  white-light  signal  processing. 

WHITE-LIGHT  OPTICAL  SIGNAL  PROCESSING 

We  shall  now  describe  an  optical  signal  processing  technique  that  can  be 
carried  out  by  a  white-light  source,  as  illustrated  in  Fig.  1.  The  white-light 
signal  processing  system  is  similar  to  that  of  a  coherent  system,  except  for  the  use 
of  a  white-light  source,  source  encoding  mask,  signal  sampling  grating,  multispectral 
filters  and  achromatic  transform  Jcnscs.  For  example,  if  we  place  a  signal  trans¬ 
parency  s(x,y)  in  contact  with  a  sampling  phase  grating,  the  complex  light  field  for 
every  wavelength  A  behind  the  achromatic  transform  lens  Lj  would  be 

E(p,q;A)  =  // s(x,y)  exp(ipox)exp[-i (px+qy) ]dxdy  -  S(p-pQ,q)  (1) 

where  the  integral  is  over  the  spatial  domain  of  the  input  plane  Pj_,  (p,q)  denotes 
the  angular  spatial  frequency  coordinate  system,  pD  is  the  angular  spatial  frequency 
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of  I  In*  iiiunpl  tug  |  >  I  •  t  f  • » •  p.i  in  ill)',,  iiml  !!(p,i|)  In  Mu*  I  M|n-i  I  Mini  ill  ii(:',y).  II  we 

write*  M«|.  (I)  In  t  lit*  liirm  nl  I  I  lion  r  npallnl  coord  lunl.o  nyulim  (it, ft),  wo  have, 


K(a.R;A>  «  SOi  -  2);  P„,tf) 


(2) 


where  p  =  (2it/Xf)a,  q  =  (2x/X£)G,  and  f  is  the  focal  length  of  the  achromatic  trans¬ 
form  lens.  Thus,  we  see  that  the  Fourier  spectra  would  disperse  into  rainbow  color 
along  the  a  axis,  and  each  Fourier  spectrum  for  a  given  wavelength  X  is  centered  at 
o  «  ±(Xf/2x)p0. 


In  signal  filtering,  we  assume  that  a  sequence  of  complex  spatial  filters  for 
various  Xn  are  available,  i.c. ,  ll(pn,q„) ,  where  pn  *  (2ii/Xnf)u,  n,.  =  (2it/Xnf )  ft.  In 
practice,  all  the  processing  signals  are  spatial  frequency  limited;  the  spatial 
bandwidth  of  each  spectral  band  f Liter  ll(pn,qn)  la  also  hand! imiled ,  such  os 


•  ,  f11  <!*„ •>!„>. <  “  <  «2 

HCp„,qn)  -  ). 

V..  0,  otherwise 


(3) 


where  =  (Xnf/2x)(p0  +  Ap)  and  a£  =  (Xnf /2tt)  (p  -  Ap)  are  the  upper  and  the  lower 
spatial  limits  of  H(pn,qn),  and  Ap  is  the  spatia?  bandwidth  of  the  input  signal 

*(x,y). 


.The  limiting  wavelengths  of  each  H(pn,qn)  can  be  written  as 


PG  +  ^ 


"  AP 


t  a>  Xt  e  A  .  » 

n  PQ  -  Ap  h  n  po  +  Ap 


The  spectral  bandwidth  of  H(pn,qn>  is  therefore 


AX  -  X  — 
n  n2 


4pQAp 


^A£ 


P*  -  (Ap)2  po  n 


(4) 


(5) 


If  we  place  this  set  of  spectral  band  filters  side-by-side  positioned  over  the 
smeared  Fourier  spectra,  then  the  intensity  distribution  of  the  output  light  field 
can  be  shown  as 

N 

I (x,y)  =  I  AXn  |s(x,y;Xn)  *  h(x,y;Xn>|  2  (6) 

n=l 


where  h(x,y;X)  is  the  spatial  impulse  response  oT  H(pn,on)  and  *  denotes  the  con¬ 
volution  oper/il  Ion.  Tlnm,  the  proponed  wli  I  l.o- 1  IghL  nlgnal  proconjior  In  capable  of 
processing  the  signal  In  complex  amplitude.  Since  the  output  intensity  is  the  sum 
•f  the  mutually  incoherent  narrow  band  spectral  irradiances,  the  annoying  coherent 
artifact  can  be  eliminated.  Furthermore,  the  white-light  source  contains  all  the 
color  wavelengths;  the  processor  is  very  suitable  for  color  signal  processing. 


55 


Sl’KCTKAI.  HAND  i-TKRINC.  ,  SH'.NAl.  SAMI’ LI  NT.  AND  SOimCK  KNCOIMNd 


As  stated  earl ler,  in  white-light  signal  processing  wc  would  approach  the 
problem  from  a  different  standpoint.  For  example,  if  signal  filtering  is  two- 
dimensional  (e.g.,  2-1)  eor  re  Ini:  Jon  operation),  wo  would  synthesize  n  set  of  narrow 
spectral  band  niters  for  each  >n  for  the  entire  smeared  Fourier  spectra,  as 
illustrated  in  Fig.  2(a).  On  the  other  hand,  if  the  signal  filtering  is  one- 
dimensional  (c.g.,  deblurring  due  to  linear  motion),  a  broadband  fan-shape  spatial 
filter,  to  accommodate  the  scale  variation  due  to  wavelength,  can  be  utilized  as 
illustrated  in  Fig.  2(b).  Since  the  filtering  is  taking  place  with  the  entire 
spectral  band  of  the  light  source,  the  artifact  noise  can  be  suppressed  and  the 
white-light  processing  technique  is  also  suitable  for  color  image  processing. 

There  is,  however,  a  temporal  coherence  requirement  imposed  upon  the  signal 
filtering  in  Fourier  plane.  Since  the  sc.nlc  of  the  Fourier  spectrum  varies- with 
wavelength,  a  temporal  coherence  requirement  should  be  imposed  on  each  spatial 
filter  at  the  Fourier  plane.  Thus,  the  spectral  spread  over  each  filter  11  (pn , <!n)  is 
imposed  by  the  temporal  coherence  requirement,  i.e., 

AX  /  A 

^.^£<<1  (7) 

11  r° 

From  this  inequality,  a  high  degree  of  temporal  coherence  is  achievable  in  the 
Fourier  plane  by  simply  increasing  the  spatial  frequency  of  the  sampling  grating. 
Needless  to  say,  the  same  tempoal  coherence  requirements  of  Eq.  (7)  can  also  be 
applied  for  a  broadband  fan-shape  filter. 

There  is  also  a  spatial  coherence  requirement  imposed  at  the  input  plane  of 
the  white-light  signal  processor.  The  spatial  coherence  function  at  the  input 
plane  can  be  shown  (refs.  14,15), 


I'(x-x')  =  //y(xo)cxpli2TT  (x-x’)]dxo  (8) 

which  essentially  is  the  Van  Cittert-Zernike  Theorem  (refs.  16,17),  where  y  (xD) 
denotes  the  intensity  distribution  of  the  source  encoding  function. 

From  the  above  equation,  We  see  that  the  spatial  coherence  and  source  encoding 
functions  form  a  Fourier  transform  pair,  i.e., 

'  Y(xo)  -^tr(x-x')]  (9) 

and 

r(x-x')  =3F-'iy(x0)]  (10) 

where,-?  denotes  the  Fourier  transformation.  This  Fourier  transform  pair  implies 
that  if  a  spatial  coherence  function  is  given  then  the  source  encoding  function  can 
be  determined  with  the  Fourier  transformation  and  vice  versa.  We  note  that  source 
encoding  function  can  consist  of  apertures  of  any  shape  or  complicated  gray  scale 
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I  real  quantity  wliieli  Is  restricted  by  tbc  following  physical  realizable  condition: 


! 


—  -0  <  Y(’xo)  <  1  (11) 

In  white-light  signal  processing,  we  would  search  for  a  reduced  spatial 
coherence  requirement  for  the  processing  operation.  With  reference  to  this  reduced 
spatial  coherence  function,  a  source  encoding  function  that  satisfies  the  physical  ' 


realizability  condition  can  be  obtained.  One  of  the  basic  objectives  of  the  source 
encoding  is  to  alleviate  the  inability  of  a  physical  white-light  source.  Further¬ 
more  the  source  encoding  also  Improves  Lite  utilization  of  the  light  power  such  that 
the  optical  processing  can  be  carried  out  by  an  extended  source. 

We  shall  now  illustrate  an  application  of  the.  source  encoding,  signal  sampling, 
and  filtering  for  a  white-light  signal  processing,  hot  us  now  consider  a  poly¬ 
chromatic  Image  subtraction  (ref.  18).  The  image  subtraction  of  Lee  (ref.  19)  that 
we  would  consider  is  essentially  a  one-dimcnsiohal  processing  operation,  in  which 
a  1-D  fan-shape  diffraction  grating  should  be  utilized,  as  Illustrated  in  Fig.  3. 

We  note  that  the  fan-shape  grating  (i.e.,  filter)  is  imposed  by  the  temporal 
coherence  condition  of  Eq.  (7).  Since  the  image  subtraction  is  a  point-pair 
processing  operation,  a  strictly  broad  spatial  coherence  function  at  the  input  plane 


Is  not  required.  In  other  words,  if  one  maintains  the  spatial  coherence  between 
the  corresponding  Image  points  to  be  subtracted  at  the  input  plane,  then  the  sub¬ 
traction  operation  can  be  carried  out  at  the  output  image  plane.  Thus  instead  of 
using  a  strictly  broad  spatial  coherence  function,  a  reduced  spatial  coherence 
function  may  be  utilized,  such  as 

T(y-y')  *  6(y-y'-ho)  +  6(y-y’+ho)  (12) 

where  2h0  Is  the  main  separation  between  the  two  input  color  transparencies.  The 
source  encoding  function  can  therefore  be  evaluated  by  through  the  Fourier  transform 
of  Eq.  (9),  such  as 

2irh 

Y(yQ)  =  2  cos(-~-  yQ)  (13) 

Unfortunately  Eq.  (13)  is  a  bipolar  function  which  is  not  physically  realizable. 

To  ensure  a  physically  realizable  source  encoding  function,  we  let  a  reduced 
spatial  coherence  function  with  the  point-pair  coherence  requirement  be  (ref.  20) 

sin(jp  1  y— y  *  | ) 

r(|y-y'|)  =  - - - Bine  (-J2L  |y-y'|>  (H) 

Nh.Ii»(-|  -  |  y-y  * } )  «» 

o 

where  N  >>  1,  a  positive  integer,  and  w  <<  d.  Equation  (14)  represents  a  seguence 
of  narrow  pulses  which  occur  at  every  |y-y'|  =  nhQ,  where  n  Is  a  positive  integer, 
and  their  peak  values  arc  weighted  by  a  broader  sxnc  factor,  as  shown  in  Fig.  4(a). 
Thus,  a  high  degree  of  spatial  coherence  can  be  achieved  at  every  point-pair  between 
the  two  input  color  transparencies.  By  taking  the  Fourier  transformation  of  the 
reduced  spatial  coherence  function  of  Eq .  (14),  the  corresponding  source  encoding 
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where  w  is  the  slit  width,  d  ■  (lf/hQ)  is  tlic  separation  between  the  sifts,  and  N 
is  the  number  of  the  slits.  Since  y ( | y | )  is  a  positive  real  function  which 
satisfies  the  constraint  of  Eq.  (11),  the  proposed  source  encoding  function  of 
Eq.  (15)  is  physically  realizable. 

In  view  of  Eq.  (15),  we  also  note  that  the  separation  of  slit  d  is  linearly 
proportional  to  the  The  source  encoding  is  a  fan-shape  type  function,  as  shown 
in  Fig.  4(b).  To  obtain  lines  of  rainbow  color  spectral  light  sources  for  the 
signal  processing,  we  would  utilize  a  linear  extended  wIiiLc-1  Iglit  source  with  a 
dispersive  phase  grating,  as  Illustrated  in  Fig.  3.  Thus  with  the  described  broad¬ 
band  source  encoding  mask,  sampling  grating,  and  fan-shape  sinusoidal  grating, 
a  color  substracted  image  can  be  seen  at  the  output  image  plane. 

RECENT  ADVANCES  IN  WHITE-LIGHT  PROCESSING 

It  would  occupy  lengthy  pages  to  describe  most  advances  in  white-light  signal 
processing.  We  would  however  have  to  restrict  our  discussion  to  a  few  recent 
results  that  are  considered  interesting.  Since  the  white-light  signal  processor  is 
particularly  suitable  for  color  signal  processing,  we  shall  provide  the  results 
mostly  in  color  images. 

We  shall  first  demonstrate  a  color  image  deblurring  result  due  to  linear  mo¬ 
tion.  Since  linear  motion  is  a  1-D  processing  operation  and  its  deblurring  filter 
is  a  point-by-point  filtering,  fan-shape  deblurring  filter  can  be  utilized 
(ref.  21).  Figure  5(a)  shows  a  color  picture  of  a  blurred  image  due  to  linear 
motion  of  an  F-16  fighter  plane.  The  body  of  this  fighter  plane  is  painted  in  navy 
blue-and-white  colors,  the  wings  are  mostly  painted  in  red,  the  tail  is  also  navy 
bl ue-and-whlte,  and  the  ground  terrain  is  generally  bluish-green  color.  From  this 
figure,  we  see  that  the  plane  is  severely  blurred  due  to  motion.  Figure  5(b)  shows 
the  color  image  deblurring  result  obtained  with  the  white-light  signal  processing 
technique.  From  this  result,  the  letters  and  overall  shape  of  the  entire  airplane 
are  more  distinctive  than  the  blurred  one.  The  river,  the  highways,  and  the 
forestry  of  the  ground  terrain  are  far  more  visible.  We  note  that  the  color 
reproduction  of  the  deblurred  image  is  spectacularly  faithful  and  coherent  artifact 
noise  is  virtually  nonexistent. 

Let  us  now  provide  a  color  image  subtraction  utilized  by  the  source  encoding 
technique  with  extended  incoherent  sources  (ref.  18)  as  described  in  previous 
sections.  Figure  6(a)  nnd  6(b)  show  (wo  color  imago  trmwipareno l«*?j  of  a  parking  lot 
nn  input  color  objects.  Figure  6(c)  shows  the  color  subtracted  image  obtained  by 
the  source  encoding  technique  with  extended  incoherent  source.  In  this  figure,  the 
profile  of  a  (red)  subcompact  car  can  be  seen  at  the  output  image  plane.  The 
shadow  and  the  parking  line  (in  yellow  color)  can  also  be  readily  Identified. 
However,  we  note  that  this  color  image  subtraction  result  is  obtained  by  two 
narrow  spectral  band  extended  incoherent  sources.  Extension  toward  the  entire 
spectral  band  of  a  white-light  source  is  currently  under  investigation. 
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.  Wc  tili.-il  I  now  1 1  I u:it  rate  our  experimental  renult  l.hal  a  miiHIcnlnr  sound 
•pert  rngrnm  (ref.  7?)  ran  t>«4  general  «**l  l>y  a  wh  I  to-  I  Iglil  |»  ro<- I  nj*  (crlmlaue,  ns 
11 1  us  l  rut  oil  In  I’lg.  7.  Figure  8(n)  shows  a  typical  1-D  intensity  modulated  speech 
signal  obtained  by  focusing  its  CRT  scanner  onto  a  moving  photographic  film  by  n 
Him  transport.  If  the  recorded  format  transparency  is  transporting  over  a 
rectangular  optical  window  at  the  input  plane  of  Fig.  7,  a  frequency  color  encoded 
sound  spectrogram  can  be  recorded  at  the  output  plane.  Figure  8(b)  shows  a  typical 
frequency  color  coded  speech  spectrogram  obtained  with  this  technique.  The  frequency 
content  is  encoded  from  red  for  high  frequency,  green  for  intermediate  frequency,  to 
blue  for  low  frequency.  This  color  encoded  speech  spectrogram  represents  a  sequence 
of  English  words  spoken  by  a  male  voice.  These  words  are  "testing,  one,  two,  three, 
four."  From  this  color  encoded  speech  spectrogram,  wc  see  that  excellent  character¬ 
ization  of  format,  variation  can  readily  be  seen.  because  of  the  use  of  a  white-light 
source,  the  artifact  noise  is  avoided.  As  compared  with  the  electronic  and ’digital 
counterparts,  its  white-light  signal  processing  technique  simplified  the  processing 
technology  and  the  system  is  rather  versatile  to  opt  rate.  Although  the  result 
provided  is  rather  preliminary,  it  is  the  first  color-coded  speech  spectrogram 
being  generated. 

We  shall  now  illustrate  a  computer  controlled  white-light  density  pseudocolor 
encoder,  as  proposed  in  Fig.  9.  The  spatial  encoding  is  made  by  multiplexing  a 
positive,  a  negative,  and  a  product  image  onto  a  black-and-white  photographic  film, 
as  illustrated  in  Fig.  10.  Figure  11  shows  a  sketch  of  the  normalized  transmittance 
as  a  function  of  gray  scale.  If  the  encoded  transparency  is  inserted  at  the  input 
plane  of  a  white-light  processor  of  Fig.  9,  then  a  density  color  coded  image  can  be 
obtained  by  color  filtering  at  the  Fourier  plane. 

Wc  stress  that  this  white-light  pseudocolor  encoder  offers  several  advantages 
over  the  digital  counterpart.  The  encoder  is  far  less  expensive  and  in  principle 
the  technique  offers  a  higher  image  resolution. 

Figure  12  shows  a  set  color  coded  image  of  a  woman’s  pelvis.  The  x-ray  was 
taken  following  n  surgical  procedure.  A  section  of  the  bone  between  the  sacroiliac 
joint  and  spinal  column  has  been  removed.  Tn  Fig.  32(a),  the  positive  image  is 
encoded  in  red,  the  negative  image  is  encoded  in  blue,  and  the  product  image  is 
encoded  in  green,  by  comparing  the  pseudocolor  coded  image  with  the  original  black- 
and-white  x-ray  picture,  it  appears  that  the  soft  tissues  can  be  better  differen¬ 
tiated  by  the  color  images  as  demonstrated  by  the  fact  that  the  image  contrast  in 
the  region  containing  the  gastrointestinal  tracts  is  evidently  superior  in  the  color 
image.  On  the  other  hand,  there  seems  to  be  a  degradation  in  the  resolution  in  the 
color  image  along  edges  of  the  hard  tissues.  This  is  perhaps  caused  by  two  reasons: 
Firstly,  high  frequency  information  may  be  eliminated  due  to  the  low  spatial 
frequency  encoding  gratings  (40  lJncs/mm  and  26.7  lines/mm)  employed.  Secondly, 
the  image  may  he  .smeared  due  to  the  film  development  process.  These  two  problems 
can  he  easily  corrected  by  selecting  higher  frequency  encoding  gratings  and  by 
gaining  more  experience  in  film  proems  Ing. 

Another  point  worthwhile  to  note  is  that  a  reversal  of  the  color  encoding  can 
be  easily  implemented  as  shown  in  Fig.  12(b),  where  the  positive  and  negative  images 
are  encoded  in  blue  and  red  while  the  product  image  remains  In  green.  This  color 
mixture  capability  could  be  beneficial  because  an  Image  in  different  color 
combination  may  reveal  subtle  features  which  are  otherwise  undetected.  For  instance, 
the  air  pockets  in  the  colon  of  the  patient  can  be  identified  more  easily  with 
Fig.  12(b)  than  with  Fig.  12(a).  Moreover,  a  wide  variety  of  other  pseudocolor 
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encoded  linages  can  a  I  ho  |>e  obtained  by  simply  nit  or  nail  hr  tin:  color  filter*  In  the 
Fourier  plane  of  the  white-light  proeennor. 

CONCLUSION 


In  conclusion  wo  would  stress  that  source  encoding  is  to  provide  an 
appropriate,  reduced  spatial  coherence  function  at  the  input  plane  so  £hat  the 
signal  processing  can  be  carried  out  by  an  extended  white-light  source,  'ihe  effect 
of  the  signal  sampling  is  to  achieve  a  higher  temporal  coherence  at  the  Fourier 
plane  so  that  the  signal  can  be  processed  in  complex  amplitude.  If  the  filtering 
operation  is  two-dimensional,  a  mul  t  l-spectral-band  2-1)  f  11.t  its  should  be  utilized. 
If  the  filtering  operation  is  one-dimensional ,  a  fan-shape  filter  can  be  used.  In 
sunmary,  the  white-light  signal  processor  is  capable  of  processing  the  signal  in 
complex  amplitude  n»  a  coherent  prneenuor  and,  on  the  oLhcr  hand,  it  suppresses  the 
artifact  noise  as  an  incoherent  processor.  Since  the  white-light  source  contains 
all  the  visible  wavelengths,  it  has  been  shown  very  suitable  for  color  signal 
processing. 
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Figure  1.  A  white-light  optical  signal  processor 


Figure  2(a).  A  multi  spectral-band  filter. 


Figures  6(a)  and  (b).  Black-and-white  pictures  of  the  input  color  objects. 


Figure  6(c),  A  black-and-white  subtracted  color  image. 
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Figure  8(a).  An  intensity  modulated  format 


Figure  8(b).  A  black-and-white  frequency  color  coded  speech  spectrogram 
This  spectrogram  represents  a  sequence  of  English  words 
"testing,  one,  two,  three,  four"  spoken  by  a  male  voice. 
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Figure  10.  Spatial  encoding. 


Figure  11.  Normalized  transmittance  as  a  function  of  gray  scale 
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